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A  study  of  molecular  motions  in  various  phases  of  carbon  monoxide 
has  been  performed  through  nuclear  magnetic  resonance  (NMR)  spin-lattice 
relaxation  time  T,  measurements  in  the  liquid  phase,  T,  and  quadrupole 
perturbed  NMR  spectra  in  the  3-phase  and  nuclear  quadrupole  resonance 
(NQR)  linewidth  and  T,  measurements  in  the  a-phase  in  a  sample  enriched 
24.5%  with  017  . 

The  tumbling  rotational  motion  in  the  liquid  phase  and  its  temper- 
ature dependence  are  obtained.  An  activation  enthalpy  of  198  ±3  cal/mole 
is  derived  from  the  T-,  measurements.  Information  on  features  of  the 
precessions!  motion  about  the  crystal  C  axis  in  the  8-phase  is  also 
obtained.  The  effective  quadrupole  coupling  constants  are  compared  with 
the  static  quadrupole  coupling  constant  obtained  in  the  NQR  measurement 
of  the  a-phase.  The  former  is  shown  to  be  three  orders  of  magnitude 
smaller  by  precessional  motion  averaging.  The  temperature  dependence 
of  the  effective  quadrupole  constant  is  also  derived.  It  is  shown  that, 
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as  the  temperature  increases,  the  average  deviation  of  the  processional 
angles  from  the  magic  angle  increase  from  3.23  t0.06°  at  61.9  K  to 
3.36  ±0.06°  near  the  triple  point.  T.  measurements  in  the  solid  3 -phase 
were  also  obtained. 

In  the  a-phase  the  quadripole  interaction  is  sufficiently  large  enough 
to  perform  a  pure  quadrupole  resonance  experiment.  The  measurements  in  the 
a-phase  provide  the  magnitude  of  the  intramolecular  quadrupole  interaction. 
Furthermore,  the  NQR  0i7  iinewidths  exhibit  a  peculiar  and  unexpected  be- 
havior with  a  maximum  around  37  K  and  a  divergence  for  T  <  28  K.  The 
effective  T  measurements  in  the  temperature  range  from  49  *   55  K  support 
the  assumption  that  the  linewidth  is  determined  by  the  molecular  motions  in 
the  a-phase.  This  behavior  has  been  analyzed  on  the  basis  of  a  model  of 
molecular  recrientational  processes.  These  processes  could  account  for  the 
residual  entropy  previously  observed  in  the  crystal.  Comparison  of  the  data 
with  the  theoretical  picture  allows  one  to  derive  information  on  the  details 
of  the  molecular  reorientations  and  their  temperature  dependence. 
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CHAPTER  I 

INTRODUCTION 

1 .  1  Physical  Properties  of  CO 

Carbon  monoxide  is  one  of  the  most  interesting  diatomic:  molecular 
crystals  and  has  been  studied  extensively  [1-121.  Besides  its  structural 
simplicity  and  its  similarity  to  other  diatomic  molecular  crystals 
[6,9,13,14],  CO  exhibits  some  peculiar  behavior  which  is  probably  related 
to  the  electric  dipole-dipole  interaction  of  the  molecules  [5,15].  Since 
the  observation  of  residual  entropy  [4],  interest  has  been  focused  on  the 
possibility  of  an  ordering  transition.  Gill  and  Morrison  extended  the 
studies  to  lower  temperatures  and  found  no  transition  occurring  [10]. 
However,  an  anomaly  in  heat,  capacity  near  18  K  has  been  reported  [16]. 
A  theoretical  calculation  of  the  rate  of  ordering  based  on  the  model  of  a 
twofold  barrier  for  end-over-end  rotation  in  solid  a-C0  found  the  rate  to 
be  unobservably  slow  [17]. 
1.  1.  1  Crystal  Structure  of  CO 

There  are  two  solid  phases,  a  and  - .  of  CO.   At  low  temperatures, 
X-ray  diffraction  studies  showed  that  the  crystal  structure  of  the 
a-phase  is  primitive  cubic  with  a  basis  of  four  molecules  per  unit  cell 
and  that  the  molecules  are  aligned  in  the  body  diagonal  directions  as 
shown  in  Figure  1-1.  However,  two  space  groups  Pa3  (T?)  [IS]  and 
P2.j3  (T  )  [1]  are  compatible  with  this  arrangement.  The  P2,3  space  group 
differs  from  Pa3  in  that  each  molecular  center  is  displaced  by  a  small 
amount  6'  along  the  direction  of  the  molecular  axis.  The  problem  concern- 
ing which  is  the  correct  structure  for  a-C0  is  still  in  dispute  [19]. 
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Figure  1-1.     Crystal   structure  of  a-CO.     The  molecules  are  aligned 
parallel    to  the  body  diagonals. 


In  the  3-phase  which  occurs  above  61.6  K,  the  crystal  structure 
determined  by  Vegard  [2]  was  that  of  B-N2,  and  subsequent  X-ray  diffrac- 
tion studies  of  b-N2  [20]  deduced  this  space  group  as  P63/mmc(D!ih).  In 
this  regard  one  considers  the  crystal  structure  of  3-CO,  the  same  struc- 
ture as  S-N2,  to  be  hexagonal  close  packed.  It  was  found  that  X-ray  form 
factors  were  equally  consistent  with  a  dynamic  model  where  the  molecules 
precess  about  the  C-axis  at  an  angle  e'  as  shown  in  Figure  1-2,  or  with 
a  static  disordered  model  in  which  the  molecules  are  randomly  distributed 
among  the  24  general  positions  of  space  group  P63/mmc  [20]. 
1.1.2  Phase  Diagram  of  CO 

The  phase  diagram  of  CO  has  been  determined  by  Fukushima,  Gibson, 
arid  Scott  up  to  a  pressure  of  1.75  Kbar  using  the  change  in  nuclear  spin- 
spin  relaxation  time  T"2  of  C13  as  an  indicator  of  the  transition  [21]. 
The  phase  diagram  is  shown  in  Figure  1-3.  The  a, 6  phase  transition  is  at 
61.6  K  at  equilibrium  vapor  pressure,  and  the  triple  point  is  at  68.15  K. 

Applying  the  Clausius-Clapeyron  equation,  —■  =  —■  ,  to  the  phase 
transition  boundaries  the  volume  increments  at  zero  pressure  are  deduced 
to  be  AVm  =2.50  cm3/mole  &nd   AVQg  =0.92  cm3/mole  [21].  While  the  change 
of  the  volume  in  N2  is  av^  =  2.5  cm3/rnole  [22],  which  is  comparable  to  CO. 


However,  iV  ,  in  •'.',,  has  a 


The  compression  experiments  found  no  new  phase  at  higher  pressures  [25], 
This  is  unlike  N2  where  a  new  phase  was  found  at  higher  pressures  with  a 
tetragonal  structure  [26]. 
1.1.3  Specific  Heat  and  Other  Thermodynamic  Properties  of  CO 

The  heat  capacity  of  CO  has  been  measured  [3,4,10]  and  is  shown  in 
Figure  1-4.  Tfie  analysis  of  the  results  is  compared  with  the  entropy  cal- 
culated from  the  band  spectrum  of  CO  in  the  gas  phase.  The  discrepancy 
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Figure  1-2.  Crystal  structure  of  3-CQ.  Open  circles  representing  the 
cage  structure  denote  hep  positions  of  molecules.  Only  the 
center  molecule  is  illustrated  which  processes  around  the 
crystal  C  axis  with  angle  o'  near  the  magic  angle. 
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Fiyure  1-3.  Phase  diagram  of  CO  in  the  P-T  plane.  The  ct-6  transition 
temperature  at  equilibrium  vapor  pressure  is  61.6  K.  The 
triple  point  is  at  63.15  K.  {Data  are  from  Ref.  [21]). 
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was  found  to  bo  approximately  equal  to  R  In  2.  This  was  described  as 
frozen  in  end-for-end  disorder  [4].  Mslhuish  and  Scott  [15]  calculated 
the  energy  for  the  oriented  and  disoriented  lattice  and  estimated  the 
Curie  temperature  for  the  order-disorder  transition.  They  found  Tc  to  be 
^5  K  for  CO.  However,  heat  capacity  measurements  on  solid  a-CO  have  been 
carried  to  2.5  K  by  Gill  and  Morrison  with  no  indication  of  the  occurrence 
of  the  transition  [10].  The  adiabatic  compressibility  xs  may  be  calculated 
from  the  sound  velocities  and  density  using  the  equation 

Xs  ■  Cp  (V*  -  f^)]-1  (!_,) 

where  V£  and  Vt  are  the  longitudinal  and  transverse  sound  velocities,  and 
p  is  the  density. 

Voitovich  et  al .  [27]  measured  V£  and  Vt  and  thereby  x  was  obtained. 
While  from  the  thermodynamic  equations  one  has 

XT  =  xs  +  H2   TV/Cp  (1-2) 

and 

xT/xs  =  Cp/Cv  (1-3) 

where  xT  is  the  isothermal  compressibility,  °,   is  the  thermal  expansion, 
V  is  the  molar  volume  and  Cp  and  Cy  are  the  specific  heats  at  constant 
pressure  and  volume  respectively.  Krupsii  et  al .  measured  the  thermal 
expansion  of  a-CO  and  applied  the  above  equations  using  x  from  [27]  and 
Cp  from  [4,10]  to  give  a  complete  table  for  lattice  parameters,  molar 
volume,  volume  expansion,  specific  heat  at  constant  volume,  and  isothermal 
compressibility  [18]. 


1 ,  2   Nuclear  Magnetic  Resonance 

Consider  an  isolated  nucleus  in  a  steady  magnetic  field  fto  and  sup- 
pose that  the  nuclear  spin  number  I  is  greater  than  zero,  so  that  the 
nucleus  possesses  a  magnetic  moment.  Quantum  mechanically  the  angular 
momentum  of  the  nucleus  is  quantized  along  the  applied  field  direction 
with  the  magnetic  number  m  taking  values  of  I,  1-1,...,  -(1-1),  -(I). 
The  Hamiltonian  can  be  written  as 

*"«  -t  ■  %  (1-4) 

where  y  is  the  nuclear  magnetic  moment  and  Ho  is  an  externally  applied 
field.  The  magnetic  moment  is  related  to  the  nuclear  spin  angular  momen- 
tum by 

vtit  =  y  (1-5) 

where  y   is  the  gyrcmagnetic  ratio  of  the  nucleus,  h  equals  £-  where  h  is 

C.T\ 

Planck's  constant.  The  components  of  the  magnetic  moment  are  thus  given 
by  the  (21+1)  value  of  y,  y(I-l)/l,. . . ,  -y(I-l)/I,  -v.     The  energy  levels 
of  the  nuclear  magnet  in  the  mag?ietic  field  Ttes  are  therefore  equally 
spaced  taking  values  of  myHo/I,  where  m  -^  I,  (1-1),...,  -(1-1),  -I. 
From  the  classical  point  of  view  the  nucleus  may  be  regarded  as  a  magnet 
dipole  v   precessing  about  the  direction  of  the  applied  field  fio.  The  rate 
of  precession  is  given  by  the  well  known  Larmor  angular  frequency 

uL  =  yHo  (1-6) 

If  an  additional  small  magnetic  field  rh  is  applied  at  right  angles  to  "fio, 
the  dipole  y  will  experience  a  torque  of  y  x  %   tending  to  increase  the 
angle  e  between  y  and  3o.  If  T\i   is  made  to  rotate  about  to   as  axis,  with 


the  same  frequency  as  the  Larmor  frequency  m,  ,  then  the  angle  0  will 
steadily  increase.  For  rotation  frequencies  differing  from  u>o  the  coup- 
ling between  jj  and  Hi  will  merely  produce  small  perturbation  of  the  pre- 
cessiona";  motion  with  no  net  affect  [28].  Hence  the  application  of  a 
radio  frequency  field  Hi  at  a  frequency  u>.  will  cause  energy  to  be  ab- 
sorbed by  the  nucleus.  The  experimental  spectrum  will  be  an  infinite 
sharp  line  at  the  Larmor  frequency  <d.  .  However,  the  effect  of  the 
neighboring  nuclei  will  produce  a  local  magnetic  field  at  each  nuclear 
site.  This  will  cause  the  line  to  be  broadened.  However,  this  broadening 
is  also  effected  by  the  distances  and  relative  directions  between  the 
neighboring  nuclei  as  well  as  their  motions.  Nuclear  Magnetic  Resonance 
(NMR)  measurements  thus  provide  useful  information  about  the  molecular 
structure  and  the  molecular  dynamics  in  the  solid  and  liquid  phases. 
For  example,  in  a  single  crystal  from  the  linewidth  in  the  solid,  it  is 
possible  to  calculate  the  second  moment  M~.  Thus  one  will  have  the 
information  about  the  separations  and  directions  between  nuclei  in  a 
unit  cell.  While  for  a  powder  sample  it  is  not  possible  to  determine 
the  directional  relation  between  nuclei,  it  is  still  possible  to  find 
out  their  separations.  Moreover,  a  dramatic  change  of  M«  may  indicate 
the  onset  of  the  rotation  of  the  molecules  in  the  solid  and  further 
reduction  of  linewidth  could  be  a  clue  to  the  diffusion  of  the  molecules 
in  the  solid.  Other  examples  include  chemical  shifts  in  the  liquid  phase 
for  identification  and  analysis  of  the  substance,  Knight  shifts  for 
studying  information  about  the  conducting  electrons  near  the  Fermi 
surface,  observation  of  chemical  reactions  processes  Prom  the  growth  of 
the  signal  strength  [28],  and  finally  the  molecular  motions  through  the 
spin-lattice  relation  time  T,  measurements  as  one  will  see  later  in  this 
experimental  study  of  CO. 
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1.  3  Nuclear  Quadrupole  Resonance 


The  interaction  of  a  charge  distribution  with  an  electrical  potential 
due  to  the  external  sources  is  given  by 

E  =  [p(r)V(r)d3r  (1-7) 

where  V(r)  is  the  electrostatic  potential.  Regarding  the  nucleus  as  a 
charge  distribution,,  and  expanding  V(r)  about  tne  origin  of  the  center  of 
mass  of  the  nucleus,  one  obtains 

E  =  V(0)|P(r)d3r  +  EYaJ^pfr)d3r  +  \  \   V^jx^pCr^r     (1-8) 

where  x  ,  a  =  1,2,3,  stands  for  x,y,z.  respectively  and 

v  =  3*  lr  ■  o  v  =  -2SL.  tf  =  o 

o    3X   ■     U      '  V.6    3x  3X0  |P   U    • 

The  first  term  on  the  right  hand  side  of  (1-8)  is  the  electrostatic  energy 
of  the  nucleus  taking  the  nucleus  as  a  point  charge.  The  second  term  in- 
volves the  electric  dipole  moment  of  the  nucleus.  Since  the  wave  function 
of  the  nucleus  has  definite  parity,  p(r)  =  p(-r)   ,   the  electric  dipole 
moment  of  the  nucleus  vanishes  [29].  One  is  only  interested  in  the  third 
term,  the  electric  quadrupole  interaction  now.  By  introducing 

and  using  the  fact  that 

L    SX  3X 

a   a  a 

one  has  for  the  quadrupole  term 

In  quantum  mechanics,  one  considers  p  as  an  operator 


(ao) 


>W)&   '  *  K=p,St„„S  ^  -  V  0-10) 

and  defines  the  operator  0/°p'  as 

^P)=i^Va-V2^(0P)(^ 

v  2  (1"11} 

=  e     K=protons     (3xaKX3K  '  Sa&r¥       ■ 

From  (1-9)  one  writes  the  Hamiltonian 

jr   =  11     V     0(°P)  M-12) 

The  Wigner-Eckart  theorem  states  that  the  set  of  matrix  elements  of  an 
irreducible  tensor  operator  A£  differs  from  any  ether  irreducible  tensor 
B  only  by  a  constant  factor  [30].  The  tensor  operator  Q^op^  then  is 
related  tc  the  tensor  operator  Q(I)  as 

C  =   e  K=protons  (3x,Kx,3<  "  6aSr2) 

=  CQ(I) 

I  I,  t-  I.I 

-  C  (3  -a  -2   --2  -  5agl2)   .  (1-13) 

Since  Q^°p  can  be  written  as  combinations  of  irreducible  tensors  like 
Q(I).  Defining  eQ,  quadropole  moment  of  nucleus,  as 

e^"<1-1  leK=prLns  (3^-rK*)|l,I>      (1-14) 

and  operating  on  the  matrix  element  mm'  in  (1-13)  for  m  =  I,  m'  =•  I 
where  I  is  the  nuclear  spin,  it  follows  that 
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In  order  Tor  C  to  have  a  physical  meaning  it  follows  that  one  needs  I>  j 
to  have  a  quadrupole  interaction.  From  (1-12),  (1-13),  and  (1-15)  it 
follows  that 


o,3 


!Q  =  6TT2tnT    ,1  vae  C  !  CiBiB  +  Va}  -«.b  I2]  (1"16> 


In  the  principal  axes  of  the  field  gradient,  one  defines  eq  =  V?7  and  the 
asymmetry  parameter 


V   -  V 
.  VXX   VYY 


whers  VXX'  VYY'  and  VZZ  are  the  Pri'nciPa"!  f*eld  gradients.  One  then  has 
from  (1-16) 


"q =  4iM?ry  ^3Iz  -  I?  +  ^Tx "  *?»  0-17) 

One  can  also  write  (1-16)  in  the  laboratory  frame  as 

(1-18) 
where 

V„  =  -  V 

V±l  =  Vxz  *  ^yz 

V±2  4^zz"Vyy)±iVxy 

and  I+  =  Ix  +  il  ,  I_  =  Ix  -  il  are  the  raising  and  lowering  angular 
momentum  operators  respectively. 

The  quadrupole  interaction  is  therefore  due  to  the  interaction  of  the 
nuclear  quadrupole  moment  and  the  electric  field  gradient  produced  at  the 
nuclear  site  by  the  external  electric  sources.  Intermodular  electric 
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sources  usually  provide  only  a  small  amount  compared  with  those  of  intra- 
molecular origin.  Nuclear  quadrupole  resonance  (NQR)  therefore  provides 
a  sensitive  method  to  study  the  molecular  motions  as  well  as  the  elec- 
tronic structure  of  the  molecules.  Comparison  of  the  quadrupole  coupling 
constant  (QCC)  of  the  free  molecule  obtained  from  microwave  studies  and 
that  from  NQR  measurements  in  the  solid  phase  provides  information  about 
the  alteration  of  the  electronic  distribution  in  the  solid  state.  NQR 
can  also  be  used  to  study  crystal  structures,  phase  transitions,  and  im- 
purity effects  [31].  Nuclear  quadrupole  interactions  usually  provide  a 
yery  effective  mechanism  for  spin-lattice  relaxation  in  the  liquid  phase. 
Thus  it  permits  a  direct  determination  of  certain  correlation  times  char- 
acterizing the  molecular  dynamics.  The  study  of  the  linewidth  and  reso- 
nance frequency  also  elucidates  the  molecular  motions  in  the  solid  phase. 


CHAPTER  II 
EXPERIMENTAL  PROCEDURE 


2.  1  The  Sample 
The  predominant  naturally  occurring  isotopes  of  carbon  and  oxygen  do 

not.  possess  nuclear  moments.  Therefore  an  MMR  study  of  CO  has  to  resort 
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to  enriched  isotopes.  C   provides  one  of  the  possibilities  [21,32,33]. 

For  a  comprehensive  study  of  spin-lattice  relaxation  and  quadrupole  reso- 
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nance  in  the  low  temperature  phase  of  CO,  0'  is  the  best  choice.  0 

with  nuclear  spin -I  =  j   [34]  has  a  nuclear  quadrupole  moment  which, 
interacting  with  the  molecular  electric  environment,  provides  the  domi- 
nant mechanism  of  spin-lattice  relaxation.  From  nuclear  quadrupole  reso- 
nance experiments  which  can  only  be  performed  with  nuclear  spin  I>1 ,  the 
quadrupole  coupling  constant  of  0   in  CO  can  be  deduced.  Therefore  a 
great  deal  of  information  about  the  CO  molecules  can  be  learned. 

One  liter  of  gaseous  CO  enriched  to  24.5%  in  0  v/as  purchased  from 
Prochem  U.  S.  Service,  Inc.,  New  Jersey.  The  sample  was  claimed  to  have 
99.9  atom  %   of  C1". 

2.  2  The  Cryostat  and  Temperature  Measurements 

The  construction  of  the  cryostat  is  shown  in  Figure  2-1.  The  con- 
densed sample  was  held  in  a  Kel-F  cavity  of  0.43  irich  in  diameter  and 
0.78  inch  long.  Surrounding  the  sample  chamber  was  an  rf  coil  wound  from 
30  turns  of  equally  spaced  #26  cotton  covered  copper  wire.  The  induc- 
tance of  the  sample  coil  was  approximately  iO  yH  which  was  used  for  NMR 
measurements  operating  at  4.098  MHz.  Another  sample  coil  made  of  ^100 


14 


15 


£-  N2  both 

N2  orHebath 
vacuum  jacket 
radiation  shield  shc-ei 


-  vacuum  or  heater 
exchanger 

1 — -  brass     can 

_| — ,  platinum  thermometer 

-  heater  wire 
sample    line 
copper    block 
Kcl-F  sample  holder 
scir.pl.;      coil 

-  sample     chamber 

i'ie   gas  heater 
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Figure  2-1.     Sample  cell   and  temperature-controlled  cryostat. 
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turns  of  #35  enamel  coated  copper  wire  with  an  inductance  of  ^100  yH  was 
used  in  the  NQR  measurements  at  ^1.1  MHz.  The  sample  coil  was  connected 
to  the  spectrometer  through  the  small  stainless  steel  sample  line.  He 
gas  was  admitted  into  the  copper  block  on  which  a  platinum  resistance 
thermometer  was  mounted.  The  He  gas  served  as  a  heat  exchanger  between 
the  sample  and  the  copper  block.  The  platinum  thermometer  calibrated  by 
the  National  Bureau  of  Standards  was  used  for  two  purposes.  First,  the 
conventional  temperature  measurement  was  achieved  by  connecting  four  leads 
to  a  Leeds  and  Northrup  Mueller  bridge.  Secondly,  the  output  of  the  gal- 
vanometer was  phase  shifted  130°  then  fed  back  to  the  power  supply  which 
controlled  the  current  to  the  heater.  With  the  Mueller  bridge  the  tem- 
perature can  be  measured  accurately  to  0.01°  K;  however,  this  experiment, 
due  to  the  uncertainty  in  the  different  readings  of  the  N  and  R  mode  in 
the  Mueller  bridge  (because  of  using  feed-back  temperature  control)  a 
0.05°  K  uncertainty  should  be  assigned.  Between  the  brass  can  and  copper 
block,  He  gas  could  be  introduced  to  serve  as  a  heat  exchanger  in  case 
temperature  equilibrium  of  the  sample  with  the  liquid  K?   or  liquid  He  was 
desired.  To  achieve  a  sample  temperature  higher  than  the  surrounding 
bath,  this  space  was  evacuated  to  10"  '   torr  or  less.  A  thin  brass  sheet 
radiation  shield  held  in  place  with  small  Kel -F  p:eces  reduced  additional 
radiation  heat  loss.  Temperatures  between  50  K  and  11   K  were  achieved 
by  pumping  the  N2  bath  and  above  77  K  the  heater  was  applied.  For  tem- 
peratures below  50  K,  liquid  He  was  used  as  a  cold  bath  with  the  heater 
on  to  maintain  the  desired  temperature.  Before  the  sample  was  introduced 
to  the  sample  chamber,  this  space  had  been  evacuated  to  4  x  10"6  torr. 
To  avoid  solidified  sample  blocking  the  line  before  it  reached  to  the 
sample  chamber,  a  vacuum  jacket  was  situated  around  the  sample  line. 
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2.  3  The  Electromagnet 
The  electromagnet  used  in  the  NMR  measurements  was  a  Varian  4012-33 
rotatable  electromagnet  with  a  12  inch  diameter  pole  face  and  a  3.5  inch 
gap.  The  magnetic  field  was  regulated  by  a  fieldial  regulated  system 
which  utilized  a  Hall  effect  probe  inside  the  magnetic  field.  The  max- 
imum magnetic  field  produced  was  around  9  KG.  Taking  long  time  stability 
and  maximum  sensitivity  of  the  NMR  signals  into  consideration,  a  7.1  KG 
field  strength  seemed  to  be  the  best  choice.  This  in  turn  produced  a 
4.093  MHz  0   NMR  signal  in  liquid  CO.  The  rotation  angle  of  the  magnet 
could  be  adjusted  to  within  1°  of  accuracy.  The  inhomogeneity  of  the 
magnetic  field  over  the  sample  volume  was  determined  in  the  liquid  phase 
to  be  M).l  gauss  which  gave  a  T|  of  about  3.5  ms.  After  long  hours  of 
operation  (^  6  hrs.)  the  fluctuation  of  the  magnetic  field  was  found  to 
be  3  x  10~4/hr.  or  less. 

2.  4  The  Spectrometer 
The  pulse  spectrometer  used  in  this  experiment  is  depicted  in 
Figure  2-2.  This  set-up  was  used  for  the  NMR  experiments  as  well  as  the 
NQR  measurements  with  only  a  few  modifications  of  the  receiver  pre- 
amplifier. The  high  stability  General  Radio  1051  frequency  synthesizer 
was  the  main  rf  source.  Since  it  is  a  heterodyne  system,  the  frequency 
was  set  to  10  MHz  higher  than  the  resonance  frequency  to  utilize  a  10  MHz 
amplifier  and  a  phase  sensitive  detector.  This  signal  when  mixed  with 
the  gated  10  MHz  reference  frequency  was  amplified  to  approximately  2V 
peak  to  peak  rf,  before  it  entered  the  transmitter.  The  transmitter  is 
adequate  to  supply  high  power  pulsed  output  and  low  power  continuous 
output.  The  magnetic  field  produced  by  the  rf  pulse  was  calculated  to  be 
^30  gauss,  which  is  large  compared  with  the  linewidths  of  the  spectra  being 
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Figure  2-2.  Block  diagram  of  the  pulsed  NMR/NQR  spectrometer. 
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studied.  The  same  pulse  which  gated  the  10  MHz  frequency  controlled  and 
gated  the  final  stage  transmitter  also. 

In  a  pulse  nuclear  resonance  experiment,  the  induction  signal  is 
typically  obscured  for  several  microseconds  after  the  transmitter  pulse 
by  the  ringing  of  the  input  turned  circuit.  This  ringing  can  be  reduced 
by  a  weak  coupling  to  the  nonlinear  input  [35].  In  this  experiment,  a 
photo-FET  Q  switch  similar  to  the  one  described  by  Conradi  [36]  was  used. 
During  the  time  the  high  power  rf  pulse  was  on,  the  D,D?  cross-diode  set 
conducted  strongly.  The  DgD.  diode-set  was  also  conducting  to  protect 
the  receiver  preamplifier  and  the  photo-FET  Q  switch.  When  the  cross- 
diode  set  D~D,  is  conducting,  point  A  of  the  x  X  cable  appears  as  a  high 
impedance.  While  the  transmitter  pulse  is  off,  D-D*  ceases  to  conduct 
and  point  B  is  a  high  impedance. 

When  the  rf   level  reaches  0.5  volt,  the  cress-diodes  D,D2  turn  off. 
This  isolates  the  transmitter  from  the  L.  C.  tuned  circuit.  Unfortu- 
nately, in  practice  the  L.  C.  tuned  circuit  and  the  transmitter  are  still 
weakly  coupled  by  the  capacitor  C,  across  diodes  D-,D9.  This  coupling 
ties  the  voltage  across  the  sample  coil  to  the  slow  ringing  down  of  the 
transmitter  final  stage  [37],  so  the  recovery  time  is  prolonge:;.  The 
negative  bias  FET  was  adjusted  just  into  the  cut-off  region.  A  pulse 
from  the  pulse  generator  turned  on  a  LED  diode  which  provides  light  to 
turn  on  the  Q  switch.  The  source-drain  path  appears  as  a  low  resistance 
(^200ft)  which  provides  a  low  impedance  path  to  ground  for  the  energy  of 
the  ringing.  The  Q  value  of  the  input  circuit  hence  was  reduced  during 
the  time  the  Q  quenching  pulse  was  on,  causing  it  to  ring  down  quicker. 
Just  as  the  ringing  decreases  to  the  order  of  the  free  induction  signal 
('MnV),  the  Q  quenching  pulse  was  turned  off  to  restore  a  high  Q  value 
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in  the  receiver  preamplifier  for  detecting  the  free  induction  signal.  The 
amplified  signal  was  then  mixed  with  the  rf  frequency  (v  +  10  MHz)  split 
from  the  same  source.  The  combined  signal  was  then  phase  sensitive  de- 
tected by  mixing  with  the  10  MHz  reference  signal. 

The  output  signal  was  observed  on  an  oscilloscope  and  also  fed  to 
the  Biomation  802  transient  recorder  which  captured  the  fast  signal  for 
slower  transfer  to  the  Fabri-Tek  1072  signal  averager.  A  POP  8/E  mini- 
computer was  interfaced  with  the  Fabri-Tek  1072  to  perform  the  Fourier 
transformation  of  the  signal  from  the  time  domain  (free  induction  decay) 
into  the  frequency  domain  (spectrum).  The  result  was  then  plotted  on  an 
X-Y  recorder. 


CHAPTER  III 
LIQUID  PHASE 

3.  1  NMR  Spin-Lattice  Relaxation 
In  the  liquid  phase,  due  to  the  fast  molecular  tumbling  motions,  the 
quadrupole  interaction  is  averaged  to  zero.     The  experimental  NMR  signal 
will,  therefore,  be  a  single  line.  The  linewidths  are  mainly  the  result 
of  the  inhomogeneity  of  the  applied  magnetic  field.  The  random  motions  of 
the  molecules  also  provide  the  mechanisms  for  the  nuclear  spin  system  to 
relax  when  perturbed  from  its  equilibrium  state.  T,  measurements  thus 
can  be  used  to  study  the  molecular  motions.  In  the  liquid,  the  energy 
differences  between  the  Zeeman  levels  are  all  equal.  The  fast  spin- 
exchange  mechanism  thus  maintains  a  Boltzmann  distribution  of  the  nuclear 
populations  among  the  energy  levels  [38],  i.e.  a  spin  temperature  can  be 
defined  and  the  relaxation  to  the  lattice  temperature  is  governed  by  a 
single  exponential  decay.  In  order  to  describe  the  macroscopic  relaxation 
time  T,  through  the  microscopic  molecular  motions,  one  considers  the  time 
dependent  Hamiltonian  «,  (t)  caused  by   the  molecular  motions.  This  ;/,  (tj 
represents  a  small  perturbation  and  induces  transitions  between  eigenstates 
of  the  Zeeman  Hamiltonian.  The  correlation  function  of  h-,    (t)  is  then 
introduced  which  expresses  a  correlation  between  the  two  different  config- 
urations of  a  nuclear  environment,  at  two  different  times.  The  rate 
equations  among  the  energy  levels  are  then  considered.  This  connects  the 
relations  between  the  experimental  relaxation  of  the  total  magnetization 
and  the  correlation  time  of  the  tumbling  motions  of  molecules. 
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Consider  h^    (t)  a  small  perturbation  among  the  eigenstates 
| a  >,..., |e  >  ,  etc.  with  eigenenergies  a, 3  etc.  of  a  spin  system  S.  The 
transition  probability  per  unit  time  for  a  transition  from  state  | 3  >  to 
state  (a  >,   assuming  the  system  S  at  time  zero  is  in  state  JB  >,  time 
dependent  perturbation  theory  gives  [6] 

Wp-za  !*<»!*!  (t)|a><«|  fl]  (f)|B>e1Mt'-t)dt  +  C.  C.  (3-1) 


where  W 


o-B 


a  a 


Taking  the  average  W   over  a  statistical  ensemble,  and 


introducing  the  correlation  function  G  „  as 

ctB 


Ga3  ^t"t''  =  <  BlV1^  a  ><   a!ffl(t')l8  * 


(3-2) 


One  uses  the  properties  that  the  random  function  <  a|#,(t)|e  >  is  station- 
ary, i.e.  GaC  depends  on  t  and  V   only  through  the  difference  t  -  V   =  t. 
It  follows  thai 


Wo6  =  |,  j*  Go6  (T)  e-1W  fc  +  c.  C. 
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(3-3) 


Since  one  usually  considers  a  time  t  >>  —  ,  it  turns  out  that  the  limits 

of  the  integration  in  (3-3)  can  be  replaced  by  +  »  and  one  has 
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1 

+  '2 
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G  0   (t)  e"lwaBT  cIt 
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(3-4) 


In  most  cases  h,  (t)  can  be  expressed  as  a  combination  of  two  irre- 
ducible tensors  of  rank  k.  Let 


*-,  (t)  ■  I   A(P)  (I)  B(P}  (t) 

p-k 


(3-5) 
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where  A  is  the  spin  part  of  the  Hamiltonian  and  B  is  that  part  of  the 
spatial  coordinates  which  is  time  dependent.  It  follows  that  (3-4)  can 
be  written  as 

K«   ■  I        h\\   <  «|A(P)|6  >  |VP)(r)  e~i:W  d,  (3-6] 


aB  p--k  V 


where 


g(P)(T)  =   B(p)(t)  B*(p)(t-t) 

is  the  correlation  function  of  B^p;(t), 

There  are   various  possible  mechanisms  of  relaxation,  (1)  the  intra- 
molecular quadrupole  interaction,  (2)  the  intramolecular  spin-rotation 
interaction,  (3)  the  intramolecular  anisotropic  chemical  shift  interaction, 
(4)  the  intermodular  magnetic  dipole-dipole  interaction  and  (5)  the 
intermolecular  electric  quadrupole-dipole  interaction.  However,  an  esti- 
mate of  the  magnitudes  shows  that  the  spin-lattice  relaxation  is  almost 
completely  dominated  by  the  intramolecular  quadrupole  interaction. 

3.  2  Intramolecular  Relaxation  f'echanisms 
The  interaction  between  the  nuclear  moments  and  the  electrons  in  the 
CO  molecule  drives  the  transition  among  the  Zeeman  levels  of  017  nuclei. 
In  principle  there  are  three  relaxation  mechanisms  which  are  of  intra- 
molecular origin.  They  are  (1)  the  intramolecular  quadrupole  interaction, 
(2)  the  spin  rotation  interaction  and  (3)  the  anisotropic  chemical  shift 
interaction. 

Consider  (1)  the  intramolecular  quadrupole  interaction.  In  Appendix 
A  it  is  shown  that 

T]   125  {  +  '    l2  t3-7-' 


cA 


e20a  . 
where  — ,^  is  the  static  quadrupole  coupling  constant,  which  one  has 

determined  in  the  NQR  measurement  of  CO.  in  the  a-phase  to  be  4.20  MHz 

and  %2   1S  tne  autocorrelation  time  of  the  molecular  tumbling  motions.  If 

one  adepts  the  value  of  tr    from  Raman  anisotropic  scattering  [39]  to  be 

x2  =  3  x  10"  3  sec.  at  a  temperature  of  77  K  one  has  i  =4.17  sec"1 

while  the  experimental  T,  at  77  K  is  220  ms.  This  gives  i  =  4.54  sec-1 

1 
In  (2)  the  spin-rotation  interaction,  the  time  fluctuation  of  the 

interactions  of  the  nuclear  spins  with  the  angular  momentum  of  the  mole- 
cule provides  a  relaxation  mechanism  for  the  nuclear  spins.  Assuming  a 
Langevin  equation  for  the  angular  momentum  and  a  rotational  diffusion  of 
Debye  type,  Hubbard  obtained  [40] 


gh-2-     2IKTt 
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V   s.r. 

where  C,  and  C„  are  the  components  of  the  spin-rotation  tensor  perpendic- 
ular and  parallel  to  the  axis  directed  from  the  center  of  mass  of  the 
molecule  to  the  resonant  nucleus.  I  is  the  moment  of  inertia  about  a  line 
through  the  center  of  mass  and  perpendicular  to  the  molecular  axis.  I  is 
related  to  the  angular  momentum  of  the  molecule  3  by  iZ  =  tvJ.  t,  is  the 
correlation  time  Tor  the  angular  momentum  of  the  molecule  and  T?  is  the 
correlation  time  for  molecular  reorientation. 

Liquids  near  their  melting  points  are  often  described  by  a  reorien- 
tation process  called  the  "diffusion  limit"  in  which  the  molecules  turn 
through  a  very   small  angle  between  collisions  and  t0  and  t,  are  related 

c  J 

by  the  following  eouation  [41]: 
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For  linear  molecules  one  has  3  C2-,.  =  2  Ci  +  Cm  =  2  C2  .     Assunrinq  the 

ST"!        "  *  J 

diffusion  limit  for  the  temperature  of  interest  one  has  from  (3-8)  and 
(3-9) 

C; 


s.r. 
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If 
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(3-10) 


Using  the  value  C  ff  =  23.1  KHz  [42]  and  t2  from  Raman  experiments  [39] 
one  finds  L-  s  r  ~  3.6  x  10~5  sec"3  which  is  negligibly  small  compared 
with  (3-7). 

Consider  (3)  the  anisotropic  chemical  shift.  The  interactions  of 
the  external  D.C.  magnetic  field  with  the  electrons  around  the  resonant 
nucleus  in  a  molecule  produces  a  magnetic  field  at  the  nuclear  site.  The 
anisotropy  of  this  coupling  under  the  molecular  motions  may  introduce 
transitions  between  energy  levels.  In  the  fast  motion  regime  one  has  [43] 


h 


a.c. 


,  .0 


(3-11) 


where  y  is  the  gyromagnetic  ratio  of  nucleus  0.  H  is  the  applied  D.C. 
magnetic  field.  o„   and  ox  are  the  parallel  and  perpendicular  parts  of  the 
chemical  shift  tensor  in  the  principal  reference  frame  respectively.   t 
is  the  correlation  time  for  molecular  tumbling  motions  described  before. 

Using  the  value  cited  by  Appleman  and  Dai  ley  [44]  that  (o„  -  oj  for 
the  0  nucleus  in  CO  is  4.6  x  10"°  and  7,1  KG  for  the  external  magnetic 
field  H  one  calculates 


=  5.6  x  10"b  sec" 


a.c 


One  notes  from  the  small  value  that  the  relaxation  due  to  the  anisotropic 
chemical  shift  interaction  can  be  neglected.  By  neglecting  the  spin- 
rotation  and  the  anisotropic  chemical  shift  the  calculated  value  of 
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4.17  sec"  for  the  intramolecular  quadrupole  interaction  is  in  good  agree- 
ment  with  the  experimental  value  of  4.54  sec  . 

3.  3  Intermolecular  Relaxation  Mechanisms 


The  spin-lattice  relaxation  may  also  arise  from  the  interactions  of 
the  017  nuclei  with  the  intermolecular  fields.  There  are  two  mechanisms 
which,  in  principle,  could  induce  transitions  in  the  absence  of  a  rf 
field.  They  are  (1)  the  intermolecular  magnetic  dipole-dipole  interaction 
and  (2)  the  intermolecular  electric  quadrupole-dipole  interaction. 

Consider  first  (1)  the  magnetic  dipole-dipole  interaction.  Following 
Bloembergen,  Purcell  and  Pound  [45]  by  assuming  spherical  molecules  and  a 

Debye  model  of  diffusion,  one  would  have  a  ~  contribution  from  inter- 

'l 
molecular  magnetic  dipole-dipole  "interaction  given  by  the  following  equa- 
tion: 

fll       3,t  ^h% 

[rj  i.d.  r  it   -air-  (3'12) 

where  2a  is  the  closest  approach  of  the  neighboring  molecules,  D  is  the 
diffusion  coefficient  of  the  liquid  and  N  is  the  number  of  molecules 
per  cm3  .  Using  D  as  2.9  x  1C~5  cm?/sec  [46]  and  the  density  of  liquid 
CO  as  0.803  g/cm3,  which  gives  N  =  1.72  x  1022/cm3  and  2a  «  3.87  x  10"8cm, 
one  obtains 

|]  .  d  -  4.5  x  10"6  sec"1  . 

Now  consider  (2)  the  intermolecular  electric  quadrupole-dipole  inter- 
action. It  is  well  known  that  the  CO  molecule  possesses  a  dipole  moment 

—  1  8 

of  0.11  x  10"  esu  [5J.  The  electric  dipoles  of  neighboring  CO  molecules 
thus  produce  a  field  gradient  at  the  017  nucleus  site.  This  is  similar  to 
the  intramolecular  quadrupole  interaction  except  now  the  field  gradient 
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is  due  to  remote  electric  soirees.  Assuming  again  a  Debye  type  diffusion 
for  the  liquid  following  the  calculations  of  Bonera  and  Rigamonti  [47], 
one  has 

9tt    e2Q?M2N 

i  n  -"  to  f-2-rr-£  <3-i3) 

i .q.  h  Da 

where  a,  N  ,  and  D  have  the  same  definitions  as  discussed  in  the  magnetic 
dipole-dipole  interaction.  eQ  is  the  quadrupole  moment  of  nucleus  017, 
and  M  is  the  dipole  moment  of  the  CO  molecule.  Taking  the  value  of  Q 
measured  by  Stevenson  and  Townes  [43]  as  -0.026  x  10~2l+cm3,  one  calculates 

'1 


ft 


(y 


r  I     -  4.4  x  10  H  sec  1    . 
lV    i.q. 

From  the  above  considerations  one  immediately  concludes  that  the 
spin-lattice  relaxation  in  liquid  CO  is  completely  dominated  by  the  intra- 
molecular quadrupole  interaction.  This  is  usually  true  for  NMR  spin- 
lattice  relaxation  of  most  substances  possessing  a  nuclear  quadrupole 
moment  in  the  liquid  phase. 

3.  4  Temperature  Dependence  of  Spin-Lattice  Relaxation  Time 

Using  the  conventional  j  -  x  -  «■  pulse  sequence ,  one  measures  the 
spin-lattice  relaxation  time  T,  in  the  liquid  phase.  Figure  3-1  shows 
the  experimental  T,  values  at  different  temperatures.  Data  were  taken 
during  a  period  of  several  months.  Due  to  the  uncertainty  in  the  T, 
measurements,  a  10%  error  should  be  assigned  to  each  point.  In  view  of 
the  cluster  of  points  taken  at  different  times,  it  shows  no  impurity  is 
effectively  absorbed  in  the  sample  during  the  course  of  experiments. 

By  assuming  the  intramolecular  quadrupole  interaction  is  the  only 
contribution  to  T,  and  using  the  NQR  measurements  extrapolated  to  zero 
temperature,  one  would  have  ~&  ^  4. 20  .MHz.  From  this  one  has  the 
following  relation  for  liquid  CO 
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5.98  x  10' 


•H  l 


(3-14) 


Using  (3-14)  and  measurements  of  T, ,  the  temperature  dependence  of  x_ 
is  shown  in  Figure  3-2.  For  comparison  the  to  values  of  N0  are  shown  in 
this  picture  [49].  Assuming  an  activation  model  [50],  such  as 


RT 


(3-15) 


where  AG  is  the  difference  in  the  Gibbs  function  between  the  initial  and 
activated  states,  R  is  the  gas  constant  and  t  is  a  constant  and  is  ex- 
pected to  be  a  relatively  slowly  varying  function  of  temperature  and 
pressure  [50].  Substituting  (3-14)  into  (3-15)  one  has 


T,  =  C  e 
i 


AG 
RT 


where  C  is  some  constant.  From  the  thermodynamic  relati 


on 


AG  =  AH  +  T 


f3AG 


(3-16) 


(3-17) 


it  follows  that 


AH 


(3-18) 


[  3  In  T]   ] 

•R      — ~yt  \ 

'..  J 
AH  is  the  activation  enthalpy,  a  measure  of  the  energy  required  to  re- 
orient a  molecule  under  conditions  of  constant  pressure.  One  notes  that 
it  is  a  reasonable  approximation  to  assume  that  the  pressure  along  the 
equilibrium  vapor  pressure  curve  is  constant  [50]. 

3.  5  Discussion 
Ewing  [7]  compared  the  infrared  spectra  in  the  gas  and  liquid  phases 
of  CO  and  suggested  a  potential  barrier  of  V  =  i  V  (1  -  cos  2e)  with 
VQ  =  120  cal/mole  in  liquid  CO. 
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Figure  3-2. 
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Temperature  dependence  of  the  autocorrelation  time  for 
tumbling  rotational  motion  in  liquid  CO  and  for  precessional 
motion  in  B-CO.  (The  "o"  data  are  for  CO  while  the  "+"  data 
are  for  N?  from  Ref.  [49].  The  autocorrelation  time  in 
liquid  CO  and  N2  is  shown  as  "•"  from  experimental  Raman 
scattering  [39]). 
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Amey  [51]  usee  the  librational  model  suggested  by  Brot  and  compared 
this  with  the  far-infrared  spectrum  of  liquid  CO  and  calculated  a  barrier 
of  190  cal/mole.  These  are  compared  with  our  experimental  activation 
enthalpy  calculated  by  least  square  fit  to  be  198  ±3  cal/mole.  Consid- 
ering the  different  types  of  experiments,  it  can  be  regarded  as  satis- 
factory. In  fact,  NMR  T-j  measures  the  barrier  height  directly  and  this 
value  is  more  fundamental  than  the  other  parameters.  The  low  value  of 
the  barrier  potential  indicates  that  CO  molecules  in  the  liquid  phase 
undergo  nearly  free  rotation. 

If  one  uses  the  experimental  results  for  the  coefficients  of  diffu- 
sion in  liquid  CO  [46 J  and  assumes  a  Debye  model  of  viscous  fluids,  the 
Stokes-Einstein  equation  yields 

T  -  W-  (3-19) 

from  which  it  can  be  calculated  that  at  77  K,  x  ^   2.8  x  10~12  sec.  This 
may  be  compared  with  the  value  of  t  =  2.7  x  10~13  sec.  derived  from  the 
T^  measurements.  In  liquid  ,\'2  it  was  found  also  that  use  of  the  Stokes 
formula  for  the  reorientational  correlation  time  overestimated  the  value 
by  one  order  ox"  magnitude  [49]. 


CHAPTER  IV 
SOLID  B-PHASE 


4.  1  Quadrupole  Perturbed  N.MR  Spectra 
.  The  hexagonal  crystal  structure  of  3-CO  has  already  been  discussed 
in  1.  1.  1.  However,  the  inability  of  X-ray  diffraction  to  distinguish 
between  a  static  disorder  and  a  free  precession  model  has  led  to  some  con- 
fusion in  the  literature  [20]. 

The  experimental  data  presented  here  are  consistent  with  the  pre- 
cession model.  It  is  concluded  from  these  experimental  NMR  measurements 
that  the  effective  nuclear  quadrupole  coupling  constant  of  017  in  8 -CO  is 
reduced  by  a  factor  of  ^  10   by  motional  averaging.  This  demonstrates 
that  the  molecular  motions  in  B-CO  are  characterized  by  a  time  shorter 
than  the  reciprocal  of  the  static  quadrupole  resonance  frequency, 
*  10~Gsec,  and  that  the  reorientations  are  not  isotropic  because  a  non- 
zero coupling  constant  exists.  The  Hamiltonian  of  the  017  nuclei  in  B-CO 
can  be  written  as 


H  «  «2  +  ^  (4-1 ) 


where  h7   represents  the  Zeeman  energy  while  h      the  quadrupole  Hamiltonian. 
is  small  compared  to  #-,  and  can  be  treated  as  a  perturbation.  To  a  first 
approximation,  the  shift  in  frequency  of  state  |rn>  can  be  written  as 


%  =  f  HwftT  vZ2(3i2z-n 


in  >        (4-2) 


where  |m>  are  the  eigenstates  of  Zeeman  energy  and  m  is  the  magnetic 

quantum  number  of  the  nuclear  spin  01  / .  V„,  is  the  electric  field 

zz 
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gradient  at  the  O17  nuclear  site  in  the  laboratory  frame.  From  (4-2) 

5 

2 


knowing  I  =  ■*  for  017,  one  has 


%,  -  1  §  vzZ  <3*2  -  *>    •  <4-3> 

In  order  to  relate  the  V   in  the  laboratory  frame  to  the  static  value  of 
the  field  gradient  Vzz  in  the  molecular  frame,  one  considers  the  pre- 
cessional  motion  shov/n  in  Figure_4-1.  In  the  case  n  -  0,  the  second 
order  irreducible  tensor  transformation  enables  one  to  have 

Vzz  =  1  (3  COs2  6"  "  1J  VZZ  (4_4) 

where  V7-,  is  the  static  electric  field  gradient  in  the  molecular  frame, 
while  e'-  is  the  instantaneous  angle  between  the  molecular  axis  and  the 
external  magnetic  field  H  .  One  also  assumes  the  field  gradient  is  of 
intramolecular  origin.  If  one  applies  the  addition  theorem,  (4-4)  can 
be  written  as 

I  (3  cos2  e'  -  1)  (3  cos2  9  -  1)  V77  •       (4-5) 


The  effective  V   would  be  given  by  the  time  average  of  the  right  hand 
side  of  (4-5) 


V_z  =  \   (3  cos2  6-  -  1)  (3  cos2  8  -  1)  Vzz 

=  \jz  cos2  e-  -  1)   (3  sin2  Y  cos2  $   -  1)  Vz7  (4-6) 

It  turns  out  from  (4-3)  and   (4-6)   that 

Av     =  r4§&~(3  cos2  e-     -  1)    (3  sin2  Y  cos2  *   -  1)   (3m2   -  2|)   (4-7) 
m       n     ibu  4 


-H& 
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Figure  4-1.  Definition  of  angular  variables  used  in  the  text  to  discuss 
the  NMR  rotation  patterns  in  3--CQ.  The  molecule  is 
described  classically  as  precessing  rapidly  about  the 
crystal  C  axis.  Tl]    is  the  field  produced  by  the  rf  coil. 
H0  is  the  externally  applied  magnetic  field.  <j>  is  the 
angle  between  H0  and  the  projection  of  the  C  axis  on  the 
plane  perpendicular  to  fl-] . 


and  for  various  m  one  can  immediately  write 


Av+(;  =  ^jjjj-  (3  cos2  o-  -  1)  (3  sin2  Y  cos2  $  -  1) 
e2Qq   


Av 


2         

,     30h  (3  cos2  6'  -  1)  (3  sin2  Y  cos2  *  -  1) 


(4-8) 
(4-8') 
(4-8") 


AvT  -  -  ^tM-  (3  cos2  e'  -  1)  (3  sin'-  y   cos2  <£>  -  1 )  . 

The  quadrupole  splitting  is  of  the  order  of  KHz  while  the  static 
quadrupole  coupling  constant  is  of  the  order  of  MHz.  This  small  non-zero 
quadrupole  frequency  might  result  from  a  value  of  6'  close  to  the  magic 
angle  9  ,  where  3  ccs2e  =  1.  One  supposes  9'  is  characterized  by  a 
distribution  function  f(e),  where  e  =  9'  -  9  ,  which  is  an  even  function 
of  e  and  normalized.  It  follows  that 

1} 

3  cos2  9'  -  1  =     f(e)  [3  cos2  (e  +  e)  -  1]  dt 
I  o 

J..2L 
2 


f  - 


f(e)  [3  cos2  (e  +  e)-l  +  3  cos2  (e  -  e)  •-  1]  d< 


f(e)  e2de 


where  e  is  the  root  mean  square  of  the  angle  e,  i.e. 


magic  angle  8  .  One  defines 


and 


T  Av±3  *  Av±5 
7.  2 


?  Av±1   ±  Av±3 
2      2 


±v<2> 


*   „<» 


(4-9) 

'"  derivated  from 

(4-10) 
(4-1 0-) 
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And  combining  (4-8)     (4-9),   (4-10)  and   (4-10').   it  follows  thai 

(2)       6e2Qq       ?/->•?  ?   ,        ,x 

v         =  "RnTh       £l     (3  sin     Y  cos2  9-I) 


and 


80h 


(1)   3e2Qq     /-,  •  ?     ->  ,> 

v     ~80h   eo  ^  sin  Y  cos  ♦  "  ^ 


(4-11) 


(4-11 


From  the  experimental  perturbed  spectra  it  is  obvious  that  in  several 
cases  a  large  portion  of  the  sample  grew  as  a  single  crystal  of  3-CO. 
The  evidence  comes  from  the  rotation  patterns  of  the  separations  of  the 
satellites  as  shown  in  Figure  4-2,  in  this  case  at  64.7  K.  The  perturbed 
spectra  at  different  angles  are  also  shown  in  Figure  4-3.  One  notes  that 
the  linewidths  are  mainly  due  to  the  inhomogeneity  of  the  applied  magnetic 

r  q 

field.  Since  the  separation  between  m  ^  +^   and  m  a  +y  and  tnat  between 

5        3 
m  =  -j  and  m  =  -j  are  symmetric,  from  the  spectra  one  can  only  determine 


and 


1  2v 


2v 


(2) 


0) 


6  e20q 


40 


-  3  e2Qq  _2 
40   h   eo 


3  sin2  v  cos2  6  -  1 


3  sin2  y  cos2  9  -  1 


(4-12) 


(4-12') 


Noting  that  |2v 

ment  of  |2v^| 
0 


(2) 


,(1) 


,(i) 


2  \?.v{    !\   while  one  measures  [2\T  '[,  the  measure- 
an  be  used  as  a  supplement  to  check  the  accuracy  of  the 


|2 


measurement.  However,  the  rotation  of  the  H  field  at  an  angle 

0  3 


$  =  |- would  give  |3  sin2  y   cos2  9  -1 1  =  1 ,  while  at  an  angle  of  9  shifted 
by  j  more  would  give  J3  sin2  Y  cos2  9  -1 |  =  |3  sin2  y   -1 | .  It  is  obvi- 
ous for  symmetry  reasons  that  the  rotation  pattern  of  (4-12')  as  shown  in 
Figure  4.2  is  a  periodic  function  of  it,  having  two  local  maxima.  As  in 
Figure  4.2  it  is  at  9  -  9Q  =  32°,  where  <po   is  a  laboratory  reference  angle, 
that  gives  9  =  |  .  At  this  angle  (4-11)  and  (4-1 1')  can  be  simplified  as 
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Figure  4-2.  The  NMR  perturbed  spectra  of  8 -CO  at  different  ($ 
angles,  where  <j)0  is  a  laboratory  reference  angle. 
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In  fact  the  values  of  [3  sin2  y  cos2  6  -  1|  as  a  function  of  <j>  have  two 
local  maxima.  For  ~  <  sin2  y  <  1  the  first  (larger)  maximum  located  at 
<J>  =  0°  while  for  0  <  sin2  y  <   -  [t   is  at  <i>  *  90°  [49]. 

Comparing  with  different  single  crystals,  it  is  possible  to  determine 
which  maximum  corresponds  to  <j>  =  90°  so  that 


3  sin2  y  cos2  4>  -  lj  =  1 


(4-14) 


It  is  clear  that  one  can  also  have  an  angle  y   such  that 


3  sin2  y  cos2  <j>'  -  1 


From  (4-14)  one  determines 


sin  y  =  ±  - 


COS  <}>' 

(O 


(4-15) 


From  the  experimental  value  of  |2vv  '\   and  the  value  of  the  static 

quadrupole  coupling  constant  v/hich  one  found  in  the  NOR  measurements  in 

a-CO  one  can  determine  the  value  of  e  .  If  one  defines  the  effective 

o 

quadrupole  coupling  constant  as 

<<*!&>.   .  SO  «   1  "" 


h  'pff  -  "F  VZZ  7   <3  cos2  e'  -1) 
eft 

=  eiQa  £2 

2h  £o 


(4-16) 


one  would  have  from  (4-12),  (4-12")  and  (4-14) 
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eff 


(4-17) 


and 


.(i) 


1 .-  eiQa  > 

>0  "   h 


eff 


(4-17*) 


The  experimental  data  from  the  separations  of  the  satellites  give 


<  =-tP  >  =  7. 


04  KHz   at   T  =  64.7  K 


eff 


Comparing  this  data  with  the  static  quadrupole  coupling  constant  of 
3.83  MHz  one  finds  that 


e!Qg.>  . 


1.67  x  10 


eff 


•3  e2Qq 
h 


(4-18! 


which  is  three  orders  of  magnitude  smaller.  Figure  4-4  shows  the  effec- 
tive quadrupole  coupling  constant  for  several  temperatures  of  B-CO.  With 
increased  temperature  one  expects  e  to  increase,  causing  the  effective 
quadrupole  coupling  to  increase  with  temperature.  This  is  consistent  with 
the  experimental  results  where  e  =  3.23°  ±  0.06°  at  T  =  62  K  and 
eQ  =   3. 36  ±  0.06°  at  T  =  68  K  near  the  triple  point. 
4.  2  T]  Measurements  in  s -CO 
It  is  shown  in  Appendix  B  that  if  all  lines  in  the  perturbed  spectrum 
are  saturated  and  one  assumes  W,  =  W«  where 


Wl  -  W5 


W2  =2W5    1 

2  -  2 


(4-19) 
(4-19") 
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then  the  spin-lattice  relaxation  is  governed  by  a  single  exponential 
decay 


1   -  4  T7 

T-,  "  5W1 


(4-20) 


In  the  course  of  the  experiment  a  sequence  of  16  pulses  was  applied 
before  a  ^  -  T  -  j    Pu'se  sequence  was  employed  to  measure  the  spin- 
lattice  relaxation  time.  It  was  found  that  within  experimental  error  the 
nuclear  magnetization  recovered  as  a  single  exponential.  Furthermore  the 
normal  |-  -  t  -  5-  pulse  sequence  would  give  an  almost  identical  result 
except  for  very   small  t.  This  may  come  from  the  small  separations  between 
all  lines  and  a  £  pulse  would  suffice  to  saturate  all  the  lines.  The 
single  exponential  decay  supports  our  assumption  of  W,  =  W«.  Following 
the  same  step  as  in  liquid  CO,  and  assuming  the  relaxation  is  dominated 
by  the  intramolecular  quadrupole  interaction,  one  has 

\2 


3 
125 


e2Qc 


(4-21) 


where  t'r   now  is  the  autocorrelation  time  for  the  precessional  motion 
around  the  crystal  C  axis.  Using  — ip  from  NOR  measurements  in  the 
a-phase  one  has  the  following  relation  between  T-,  and  t': 


T2=  5. S3  x  10 


-14  1 


|4-22) 


Shown  in  Figure  3-2  is  the  temperature  dependence  of  x'0   and  T„.  Although 

C.  L. 

there  is  a  different  meaning  for  "C  and  T?,  they  specify  the  correlational 
motions  in  B-phase  and  liquid  CO  respectively.  Figure  3-2  shows  that 
there  is  a  discontinuity  at  the  transition  temperature.  It  is  reasonable 
that  the  characteristic  time  of  the  motions  in  e-CO  takes  a  longer  period 
than  that  in  liquid  CO,  while  the  discontinuity  in  the  opposite  sense  of 
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3-N?  near  the  transition  temperature  is  not  clear  at  this  moment  [52]. 
For  3-N«  the  same  assumptions  were  made  as  one  did  in  e-CO,  namely  that 
W,  =  W?  and  that  all  lines  are  saturated.  One  therefore  has 

1  =  5W,  (4-23) 

'l     ' 


in  liquid  N?  as  well  as  in  g-Np. 


CHAPTER  V 
SOLID  a- PHASE 

5.  1  Mgjecular  Properties  of  a- CO 
Solid  a-CO  is  one  of  the  simplest  diatomic  molecular  crystals  and 
has  been  an  interesting  subject  for  studying  the  lattice  dynamics  of 
molecular  crystals  [6,11,12,53-55].  The  crystal  structure  has  been  dis- 
cussed in  1.  1.  1.  One  notes  that  in  the  face  center  cubic  crystal 
structure  the  locations  of  the  molecules  and  the  different  orientations 
along  the  diagonal  can  be  related  to  the  molecular  dipole-dipole  and 
quadrupcle-quadrupole  electrostatic  energies  for  the  various  pair  orien- 
tations [15].  The  specific  heat  study  of  CO  by  Clayton  and  Giauque  [4] 
observed  a  residual  entropy  of  approximately  R  In  2   at  low  temperatures 
of  a-CO.  Gill  and  Morrison  [10]  extended  the  measurements  down  to  2.5  K 
with  no  indication  of  the  occurrence  of  an  order-disorder  transition. 
However,  it  can  be  mentioned  that  a  small  anomaly  has  been  observed  in 
the  heat  capacity  of  a-CO  near  18  K  and  ascribed  to  the  freezing-in  of 
molecular  head-to-tail  reorientation  [16]. 

5.  2  Nuclear  Quadrupole  Resonance  in  a-CO 
To  study  the  molecular  dynamics  in  a-CO,  one  can  employ  the  017 
nuclear  quadrupole  resonance,  since  the  quadrupole  coupling  constant  is 
large  enough  to  be  able  to  observe  the  NQR  frequency  directly.  For  017 
in  a-CO  intermolecular  contributions  to  the  field  gradient  are  \-ery 
small  compared  to  the  intramolecular  interactions.  Each  molecule  can 
therefore  be  regarded  as  an  isolated  spin  system.  The  electric  field 
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13 


gradient  (t'FG)  tensor  at  the  nuclear  site  is  thus  axially  symmetric  and 
the  molecules  in  a-CO  differ  little  from  that  in  the  gas  phase.  The  quad- 
rupole  Hamiltonian  can  be  written  as 

VIi^Bt)  l"31!-:2]  (5-1) 

where  eq  (T)  is  the  time  average  of  the  principal  component  of  the  EFG 
tensor  taken  over  the  vibrational  states  at  temperature  T  and  eQ  is  the 
nuclear  quadrupole  moment.  It  follows  that 

»5     3  ■  4*m  (5-2! 

±2"  -^±2  " 


=    J   e2Qg(T) 

4  «4     20  %  (5-3) 

The  zero  temperature  NQR  frequency  would  in  turn  determine  the  static 
quadrupole  coupling  constant  (QCC).  Since  by  extrapolating  the  NQR  fre- 
quency to  0  K  one  gets 

^5    3=1 .15  MHz 

±2  <-±2- 

With  zero  point  correction  estimated  from  Raman  spectra  [55],  one  obtains 


e20q 


4.20  MHz 


r=  u'^:  mz  (5-4) 

The  measurements  of  the  NQR  frequency  versus  the  temperature  are  shown  in 
Figure  5-1.  The  static  QCC  was  also  used  to  derive  t?_   and  Tx  in  the 
liquid  and  3-phases  respectively.  The  microwave  experiment  gave  4.43  MHz 
f0r  __£   This  was  refined  by  Flygare  and  Weiss  to  be  4.48  MHz  by  con- 
sidering in  addition  the  spin-rotation  interaction  in  the  microwave 
spectrum  [56].  Since  the  same  averaging  factor  due  to  molecular  stretch- 
ing motions  may  exist,  in  both  solid  and  liquid  phases,  the  experimental 
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Figure  5-1.  Temperature  dependence  of  the  pure  quadrupole  resonance 
of  017  in  a-CO. 
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NQR  determined  QCC  from  the  a-phase  should  be  more  appropriate  for  deriving 
the  correlation  time  of  the  molecular  tumbling  motions  in  the  liquid  phase. 
5.  3  NQR  Linewidt.hs  and  Molecular  Motions 

The  most  interesting  thing  that  occurs  in  the  NQR  measurements  of 
the  a-CO  is  the  anomalous  behavior  of  the  linewidths.  Figure  5-2  shows 
the  temperature  dependence  of  the  linewidths  and  some  T,  measurements, 
while  Figure  5-3  shows  a  three  dimensional  representation  of  the  spectra. 
As  one  can  see  there  is  a  maximum  at  ^  37  K  and  a  divergence  below  %  28  K. 
One  might  speculate  the  broadening  of  the  linewidths  is  related  to  static 
effects  such  as  anomalies  in  the  expansion  coefficient,  metastable  states, 
strains,  etc.  However,  the  thermal  expansion  of  a-CO  shows  no  anomalous 
behavior  in  this  region  [18]  and  the  broadening  persists  for  long  periods 
of  time  K  5  hours).  As  one  will  see  later,  the  spin-lattice  relaxation 
measurements  in  the  temperature  range  48  *  60  K  seem  to  support  the  hypo- 
thesis that  the  linewidth  is  actually  dominated  by  the  motional  contri- 
bution. 

The  occurrence  of  orientational  disorder  in  the  a-CO  crystal  has  been 
expected  since  the  first  estimate  of  the  residual  entropy  was  made  [4]. 
However,  the  frozen-in  static  disorder  cannot  explain  the  peculiar  be- 
havior of  the  linewidth  observed.   In  the  light  of  the  temperature  be- 
havior of  the  017  NQR  linewidth  and  spin- lattice  relaxation  rates,  it 
will  be  assumed  that  the  disorder  in  the  orientations  of  the  CO  dipoles 
is  of  a  dynamic  nature  at  least  for  T  £28  K. 

In  the  dynamic  model  one  assumes  that  the  molecules  execute  libra- 
tional  motions  around  the  equilibrium  positions  plus  sudden  reorientations 
among  the  various  directions  which  are,  in  principle,  permitted  by  the 
crystal  structure.  In  order  to  explain  the  peculiar  behavior  of  the 
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Figure  5-2.  Temperature  behavior  of  the  017  NQR  linewidths  in  a-CO. 
The  insertion  shows  some  effective  T^  measurements. 
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1  ine.vidth,  that  the  quadrupole  resonance  is  not  affected  (normal  Bayer 
type)  by  the  reorientation,  one  assumes  the  molecule  executes  Vibra- 
tional motion  plus  a  sudden  near  tt  reorientation  [57].  In  presence 
of  motions,  the  quadrupole  Hamiltonian  can  bo  divided  into  two  parts. 

HQ  =  <  HQ  >   +  ^Q  "  <  HQ  ^  '  ^'^ 

The  average  Hamiltonian  <  Hr>  gives  the  NQR  frequency  while  the  flue- 
s' 

tuating  part  drives  the  relaxation  process  and  causes  the  line  broad- 
ening [58].  For  017  in  the  CO  molecule,  one  assumes  a  rigid  EFG  tensor 
of  intramolecular  origin  with  cylindrical  symmetry  which,  as  a  conse- 
quence of  motions,  moves  its  principal  axes  in  the  laboratory  frame.  A 
straightforward  tensor  transformation  gives  for  the  instantaneous  EFG 
components  (using  s(t)  as  defined  in  Appendix  C) 

Vxx  (t)  =  -  3J  (1  -  3*|)  (5-6) 

Vyy  (t)  .  -  S|  [i  -  lt£  (l-s(t))  -  3*2  -  ^f-   (1  -  sit))] 

Vzz  (t)  ■  -  ^  [-2  +  3*i  +  3*?i  +  2£l  (l  -  S(t))  +  ^  (1  -  s(t))] 

Vxy  (t)  =  'H1^   C-  T(1  "  s(t))  +""  s(t)] 

vxz  (t)  =     t9-^  [V  (1  "  s(t))  +  s(t)  +~r  (1  '  s(t))] 

Vyz  (t)  =  -  ^9.  [  4J  (i  _  S(t))  +  ♦„  ] 

where  4>„  and  (j>x  are  the  in-plane  and  out-of-plane  librational  angles  with 

respect  to  the  -n   reorientation  respectively.  £*  is  a  random  variable  of 
i 

order  <  <{>£  >2  which  allows  reversal  of  an  angle  which  is  not  exactly  n  as 
discussed  in  [57].  Taking  the  librational  average  of  (5-6)  and  assuming 
no  correlation  between  <f>„,  4>x  and  <j>*  one  has 


and 
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<Vxx  >^  -S|0   -  3<$£>)  C5-7J 

<  V      >  =  -  e4  [1   -  &J£    (1   -  <  s(t)  >  )  -  3  <  tf  >  ] 

<  Vzz  >  -  eq  [1   -  |  <  tf  >  -  |  <^  >   - •  34J1  (  I   -  <  s(t)  >  )] 


<  vxy  >  -  <  vxz  >  -  <  vy2  >  -  0 


(5-8) 


In  the  casa  of  disorder,  where  <  s(t)  >  =  0,  one  has 

<  Vzz  >  =  eq  (1   -  j<$l  >  -  2  <  *"  >  "      aT  ^     • 

5        3 
The  resonance  frequency  between  m  ~  ± -*  and  m  = ± j    will  be 

3  e2Qa  M   3  x2  3    2    3  *2%         /eg* 

VQ  =  TO   h   (1  "  2  *  *S  *  "  1  K  ^   >   4  *  j   '       lb  y; 

This  is  the  usual  Bayer  type  NQR  frequency.  One  also  notes  that  for  the 
order  phase  <  s(t)  >  f   0 

vQ  -  t§  ^ga  [i  -  f  <  «2  >  -  f  <  a  >  -  -3jf  <i  -  <  set)  >)] 

(5-10) 
which  is  higher  than  in  (5-9). 

One  also  notes  that  if  reorientation  is  sufficiently  different  from 

ir,  i.e.  <j>*  is  a  large  angle,  vQ  would  be  significantly  affected.  In 

order  to  evaluate  the  linewidth  one  considers  the  fluctuating  parts  of  the 

components  of  the  EFG  tensor  from  (5-6)  and  (5-7)  and  gets 


'">2 


V  -  <  V   >  -  ^|g-  (*i  -  <  *i  >  ) 
xx     xx     2  rx  rx     ' 

V  -  <  Vyy  >  =  ^  [♦?  -  <  *?,  >  -  *£  (s(t)  -  <  s(t)  >  )] 


(5-n) 


V   -  <  V   >  = 
zz     zz 


2 


[<pi  -  <  *1  >  +  *?.  -<<{'"> 


+  C-(s(t)  -  <  s(t)  >  )] 


Vxy"  <  Vxy>  =  -^*x  C-^(l  -  s(t))  -♦.  S(t)] 

v  -  <  v  >  =  M  ^  s(t) 


xz 


xz 


V   -  <  V   >  = 

yz    yz 


3eq 


I9-  [  *2  (1  -  S(t))  +  ♦„]   . 


There  are  two  parts  contributing  to  the  linewidth  of  the  quadrupole 
resonance,  i.e.  the  adiabatic  terms  and  the  nonadiabatic  contribution  due 
to  the  lifetime  of  the  energy  levels  [58].  One  writes 


i  .i  +; 

T2  T'2       2 


1   +1 

T+5   T+3 
±2   _2 


(5-12) 


where 


1 

T2 


)'(t)t0'(t-T)  dT 


is  the  adiabatic  contribution  and  w*(t)  is  the  departure  of  the  instan- 
taneous resonance  from  its  average  value  and 


7  W  P 


(5-13) 


5    ,5  ±75-  ,  m 
tj   mfj     2 


\  -  1 


W. 


(5-14) 
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where  W    is  the  transition  probability  per  unit  "time  from  state  m  to 
state  rrr  by  the  fluctuation  terms  in  (5-11). 

By  neglecting  for  simplicity  the  terms  which  do  not  involve  the  re- 
orientations one  would  have 


9e2Qq 
20h 


1nf^  h  <°> 


(5-15) 


where  Jn  (to)  is  the  Fourier  transform  of  s(t)s(t  -  t).  Assuming  a  har- 
monic oscillator  equation  for  ^(t)  and  an  exponential  correlation  function 
s(t)  with  characteristic  time  t^',  one  has 

1     f  9e  Qg  I2  T  5  <$*2>     ty  .  20  <  <px2>    xl' 


+  1 Jn  (j  <o  ) 


15-16) 


where  xy   is  the  correlation  time  for  reorientation,  to  is  the  resonance 

5        3 
frequency  between  m  = ±  f  and  m  =  ±j.     wx  is  the  frequency  for  the  harmonic 

oscillator  $j.(t)  which  is  much  higher  than  to,  while 
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Combining  (5-15),  (5-16)  and  (5-17)  one  obtains 
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(5-18) 


1  3 

In  the  fast  motion  regime,  where  J™  (0)  ~  J™  [j  a)   -  Jj  (w)  -  J^  (tw), 

one  can  neglect  terms  involving  <  $**+>,   so  that  (5-18)  reduces  to 
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r  9e2Qq  1 
20n     _ 
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<(J)*2> 

r  s- 

.       ^ 

9 
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2 

4  +  u>*t£' 

(5-19) 


From  (5-19),  which  has  a  maximum  at  xX'  -  1.188  -  ,  one  notes  that 

at  ^37  K  the  correlation  time  for  reorientation  of  CO  is 

6  5 

x"  -   1.054  x  10  sec.  The  spin-lattice  relaxation  for  NQR  with  I  =  j  is 
2  *■ 

not  governed  by  a  single  exponential  decay  [59].  One  can  define  an 

r  i ' 


effective  relaxation  rate 
have 

'1 


eff 


which  in  light  of  (5-18)  one  would 


1 
T 


eff 


(5-20) 


T 


eff 


where  a  is  a  constant  of  the  order  of  unity.  The  measurement  of 

in  the  temperature  range  of  49 ^  60  K  supports  the  hypothesis  that  the 

linewidth  is  due  to  the  motions.  In  fact,  for  T  -  55  K  one  has 

\  =  0.25  x  103sec_1  and  according  to  (5-20)  the  linewidth  should 

'lj  eff 
be,  for  a  =  1,  0.25  KHz  in  satisfactory  agreement  with  the  experimental 

result.  One  notes  that  a  background  contribution  to  the  linewidth  can  be 

expected  from  dipolar  interactions,  strains,  and  defects,  etc.  For 

T  <  37  K  the  linewidth  seem  to  point  out  the  changeover  to  a  slow  motion 

reoime  where  Sv  decreases  on  cooling.  For x£»—  one  would  have  from  (5-19) 


9    **2    1 
2ir         Tp 


(5-21) 


It  should  be  pointed  out  that  in  the  slow  motion  regime  the  perturba- 
tion approach  implicit  in  the  above  derivation  breaks  down.  Alexander  and 
Tzalmona  [60]  showed  for  slow  motions  in  pure  quadrupole  resonance  spin- 
lattice  relaxation  time  and  therefore  the  reversal  of  the  linewidth  is 


directly  proportional  to  the  motional  correlation'  times.  It  can  be  spec- 
ulated that  the  lack  of  symmetry  of  5v  versus  T  around  T  -  37  K  is  related 
to  the  breakdown  of  the  perturbation  approach  in  the  slow  motion  regime. 

For  temperatures  lower  than  32  K,  6v  exhibits  a  dramatically  increased 
broadening  in  a  few  degrees  up  to  10  KHz  and  the  NQR  signal  is  no  longer 
detectable  below  ^28  K.  Careful  attempts  for  various  temperatures  down 
to  4.2  K  over  a  wide  range  of  frequency  searching  have  been  unsuccessful. 

A  possible  origin  for  the  behavior  of  the  linewidths  could  be  a 
change  of  nature  of  the  reorientations.  In  particular,  the  slowing  down 
of  motion  could  favor  reorientations  from  a  body  diagonal  to  an  equiv- 
alent one.  As  discussed  by  Rigamonti  and  Brookeman  [CI],  for  a  reorien- 
tation effectively  different  from  ir,  one  would  have 

I  Mv^  w?  x'-  (5-22) 

'l  £ 

where  t«'  is  the  fluctuation  time  between  the  diagonal  equilibrium 

positions. 

This  provides  one  of  the  explanations  of  not  being  able  to  observe  the 
NQR  signal  below  ^2SK.  Another  possible  origin  can  be  found  as  in  a-N„ 
where  the  broadening  of  lines  was  characteristic  of  all  mixtures  of  isotopes 
and  is  a  consequence  of  the  statistical  randomness  of  the  environments  of 
a  given  type  of  molecule  that  produces  a  dynamic  disorder  and  a  local 
inhomogeneity  in  the  average  field  gradient  [62]. 

5.  4  Discussion 

As  shown  in  Figure  5-4,  the  temperature  dependence  of  resonance  fre- 
quency of  a-CO  is  nearly  like  that  of  a-N?,  furthermore  the  reduced  NQR 
frequency  vQ  (T)  /  vQ  (0)  has  the  value  of  0.84  for  CO  as  compared  with 
0.82  for  Np  at  the  transition  temperature. [49]. 
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The  peculiar  behavior  of  the  linewidth  in  the  a-phase  NQR  measurements 
elucidates  the  possibility  of  the  motions  in  a-CO.  This  is  consistent  with 
the  observation  of  the  residual  entropy.  Some  (T^^.  measurements  also 
support  the  assumption  that  motional  effects  determine  the  linewidth. 

The  temperature  behavior  of  Sv   between  40  *  52  K  is  satisfactorily 
fitted  by  the  law 

_E_ 

and  suggests  an  activated  temperature  behavior  for  t«%  with 
E  *  1.2  K  cal/mole. 

The  assumption  of  hindered  rotation  can  be  ruled  out  since  it  would 
affect  resonance  frequency  by  a  factor  of  1  and  this  would  conflict  with 
the  microwave  measurements  of  QCC  =  4.4  MHz. 

It  could  also  be  noted  that  no  reorientation  by  an  angle  significantly 
different  froni  it  could  be  present  in  the  fast  motion  regime  where 
(T  >  37  K),  since  the  averaging  effect  would  change  vQ,  and  this  is  not 
observed  in  the  experiment.  However,  for  temperatures  T  <  32  K,  it  is 
possible  that  the  reorientation  involves  angles  different  from  rr  since  in 
the  slow  motion  regime  ut£'  >> 1  and  no  average  is  taken. 

The  dynamics  driving  the  molecular  reorientation  should  have  a  co- 
operative character  through  the  molecular  dipole  and  quadrupole  interac- 
tions. The  correlation  time  t~  introduced  above  has  to  be  considered 
a  local  correlation  time  driven  by  the  cooperative  relaxational  time  at 
a  given  wavevector.  Gill  and  Morrison  [10]  measured  the  dielectric  con- 
stant of  a-CO  down  to  6  K  and  found  no  critical  effect.  This  seems  to 
suggest  that  the  characteristic  wavevector  of  the  cooperative  excitation 
should  be  far  from  zero  and  the  possible  phase  transition,  if  any,  should 
be  of  anti- ferroelectric  type. 


APPENDIX  A 

QUADRUPOLE  RELAXATION  THROUGH  MOLECULAR 

TUMBLING  MOTIONS  IN  LIQUID  CO 

5 
In  the  liquid  phase  for  nuclear  spin  I  =  j,   the  Zeeman  energy  levels 

are  all  equally  spaced.  In  the  absence  of  an  rf  field,  transitions  intro- 
duced by  the  quadrupole  interaction  are  shown  in  Figure  A-l . 

The  actual  values  of  the  transition  probabilities  depend  on  the  de- 
tailed forms  of  the  spectral  densities  of  the  field  gradients.  However, 
all  the  upper  transition  probabilities  for  Am  =  1  are  simply  related  to 
each  other.  This  is  also  true  for  Am  =  2.  The  quadrupole  Hamiltonian 
can  be  written  as  in  (1-18) 

*Q   41(21-1)  L  3  V-Xz   l  '  '  V+  U-  lz   +  Az  l-> 

.+  V_  (I+  Iz  +  Iz  I+)  +  V+2  (IJ2  +  V_2  (I+)2  ]        (A-l) 

where  I+  and  I_  are  the  raising  and  lowering  angular  momentum  operators 
respectively.  V  represents  the  irreducible  field  gradient  tensor  com- 
ponents which,  are  related  to  the  cartesian  tensor  components  by 

V„  =  -  [3V   -  (V   +  V   +  V  )]  (A-2) 

o   yg-  L     zz   v  xx   yy    zz'J  v  ' 

V±l  ■  Vzx  *  1Vzy 

V±2  '  7  <Vxx  "  V  *  1Vxy    ' 

One  notes  that  the  V. .  are  symmetric  and  V   +  V   +  V   =0.  For  the 
raising  and  lowering  operators  one  has 
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Figure  A-l .  Quadrupole  induced  transition  probabilities  among  the 
Zeeman  energy  states  for  spin  I  -  5/2. 


I+I,  I  I,m  >  =  m  /  I(I+ll  -  m(m±l)         I.wttl  > 


and 


I-I+   I  I,m  >  =  (m±l)   /ITL+VJ  -~m[m±T7      I,m±l  >      . 

z   -     j  | 


(A-3) 

(A-3') 


One  then  calculates  the  matrix  elements  from  (A-3)  and  (A-3') 
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(A-4') 


and  also 


2'  "2 


Vz +  y±  14 


=     0 


(A-4") 


For  the  transitions  involving  Am  =±2,  one  obtains 
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(A-5-) 


From  (3-2),    (3-4),   (A-4)  and   (A-5)   it  follows 


n40  j 


>L  .  f  & 


W. 


ti40  J 


30 


SO 


e"uLT     OTTVJ,  :(t  -  t)  dt 


ei2'°LT     V+2(t)   V*2   (t  -  t)  dt 


(A-6) 
(A-61 


J-o 


where  W]   =  Wg     3  and  W2  -  2Wg     r     One  will   need  to  calculate 
V  1  V  2 

VjTtTY^    (t  -  t)     and  \^PJ*2  (t~-~7F.     To  do  so  one  considers  the 

irreducible  tensor  transformation  as  following.     For  irreducible  tensors 

k  t 

V  of  rank  k,  the  2k  +  1  components  of  V  are  transformed  accordinq  to 

P  3 

the. irreducible  representation  D.  of  the  rotation  group 


V«  =  I     vj  Dk  (a,  3,  y) 
q   t   p   pq  v  ,  p,  tv 


(A-7) 


where  (a,  B,  y)  are  the  Euler  angles  of  the  rotation  taking  the  unprimed 
reference  frame  into  the  primed  frame. 

The  Euler  angles  are  defined  by  three  consecutive  rotations:  (1)  A 
rotation  of  angle  a  about  the  z  axis,  (2)  A  rotation  of  angle  6  about  the 
y'axis  and  (3)  A  rotation  of  angle  y  about  the  z"  axis  as  shown  in 
Figure  A-2. 

The  unitary  operator  D  (a,  8,  y)  can  be  written  as 


D  (<*,  B,  y)  =  e 


Using  the  fact  that 
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Fiyure  A-2.  The  Euler  angles.  First,  a  rotation  of  angle  a  about  the 
Z  axis.  Secondly,  a  rotation  of  angle  3  about  the  Y'  axi< 
Finally,  a  rotation  of  angle  y   about  the  T"   axis. 
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which  can  be  proven  as  follows.  Let  e~  a  z  j  a  >  -'-  b  >  where  {ja>}  is 
a  complete  set  of  eigenkets  before  the  rotation  of  a   about  the  z  axis, 
{|b>}  is  another  complete  set  of   eigenkets  in  the  coordinate  system 
after  rotation.  One  has 


<  b'   J.   b  >  =  <  a"   J   a> 


=  <  a' 


a  > 


(A-10) 


It  follows  that 


V        y  e 


(A-ll) 


Using  (A-ll)  and  by  expanding,  it  can  be  proven  that 
e-iBV  =  e-iaJz  e"lBJy   e1aJz 


(A-12: 


Following  the  same  procedure,  one  has 


■iYJ 


z"  =  e" 


iaj7  _-1$Jv   -iYJz  „i6Jw   iaJ 


Jy    e 


e,MUy  e 


(A-13) 


From  (A-S),  (A-12),  and  (A-13),  it  can  be  shown  that 
D  (a,  e,  y)  =  e" 


•iaJz   -iBJy   -iYJ2 


(A-14) 


Let  |j,m  >  be  a  simultaneous  eigenket  for  J"  and  J   .   the  matrix  elements 

DJ     ,  then  will   be  written  as 
m,m 


DJ     ,  =  <  j,m 
m,m  ° 


=  <  J5m 


D   (a,   B,   Y) 


J.rrT  > 


-iaJ7        -iaJw        -iYj, 
e        z     e        Ye        z      j,m 


-lam       -iYm        ,j 
e  e  a       . 

m  m 


(A-15) 


where  the  definition  of  dji     ,  is  obvious, 
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From  the  table  [30]  one  has  for  elements  of  d2  .(g)  as 

mm 


C  =  d%   =  cos4  {%) 

22     -2-2  c 


(A-16) 


d21  =  ~d12  =  "d-2-l  =  d-l-2  =  I  Sil1  3  (i  +  C0S  &  ) 


d?0  -  d02  =  d-20  =  d0-2  =/3/§  Si"2  S 

d2-l  =  d?-2  ="d-21  =  -d-12  =  \   S1'n  B  (C°S  B  "  1} 

d2   =  d2   =  sin4  (|) 

2-2     -22  X<T 

dll  =  d-]-l  =  1  (2  C0S  6  "  1}  (C0S  3  +  1} 

di-i  =  d-n  =  ]•  (2  cos  3  +  1)  (1  -  cos  3) 

dTo  =  doi  =  "doi  =  -d?o  =  -,/57?  sin  e  cos  B 


d2Q  =  \  (3  cos2B-  1) 


In  the  principal  reference  frame  since  V..  =  0  for  i  f   j  one  would 
have  from  (A-2)  that 


0   /6  VZZ 


V..T  =  0 


V±2   2  ^VXX  "  VYY^ 
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Transformation  from  the  principal  frame  of  the  field  gradient  to  the 
laboratory  reference  frame,  yields 

v  n  =  -  V77  e"1Y  /3/2  sin  3  cos  3 
+1   /§■   ZZ 

+  1  (Vn  -  VyY)  e*i2a  e"iY  [-  |  sin  3  (1  +  cos  6)] 

+  1  (Vxx  -  vYY)  ei2a  e"iT  [-  \  sin  3  (cos  3  -  !)]   (A-17) 


V  ,  =  -  V7V  e1Y  (-  /3/2  sin  3  cos  3) 

_1   /6~  ZZ 

+  \  (Vxx  -  VYy)  e"i2a  eiY  [1  sin  3  (cos  3-1)] 

a.  1  fv          \i     ^  J2a  p1y     fy  sin  s   (co3  B  +  ^           (A-17-) 

2  ^   XX  "   VYY;  e 


V+2  =    f    VXX  e_1'2Y     W*    Sin2  ^ 
vb 


2       XX  "   VYY     e 


-l"2a      -i2Y    _t    ,3\ 
e  cosH   vTs-j 


+  1  (Vxx  -  VyY)  ei2a  e"i?Y     sin"   (f) 


(A-18) 


V-2  "    f    VZZ 
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+  £  (Vxx  -  VyYJ     e         e  sin^   (|) 


+  I(VXX-   VYY>     ^^  e1?Y     COsh{2] 
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If  one  assumes  isotropic  motions  in  the  liquid  phase  and  takes  the 
average  over  all  possible  orientations  one  obtains 
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where  n  =  — % —    is  the  asymmetry    parameter.    Assume 
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Vi   (t)  V*  (t  -  t)  =  | VTTtp^l   e"  x2 


(A-21) 


for  i  =  ±1,  ±2    where  t«  characterizes  the  correlation  function  of  V^   (t) 
Using  (A-6)  and  (A-6').   one  finds 
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In  the  fast  motion  regime  m.Tp  <<  1 
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The  relationship  between  the  experimental  T,  measurement  and  W-, ,  Wp 

are  obtained  by  considering  the  rate  equations  governing  the  populations 

among  the  energy  levels. 
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The  number  of  nuclei  in  state  m  has  the  following  time  dependence: 


dNm      I 


-a?  =  r*   ("Nm  \m-  *  V  wm:ra   '  (*-24» 

where  W  -   is  the  transition  probability  per  unit  time  from  state  m  to 
state  m't  and  N  are  the  number  of  nuclei  in  state  m.  One  can  write  down 
all  the  equations 

dN5/2  1 

"dP  =  "N5/2  Wl  "  2   N3/2  W2  +  N3/2  Wl  ^  +  A> 

+  \  N1/2  W2   (1   +  2A)  (A-25) 

dN 

■iF    =  "  I  N3/2  Wl   "N3/2W1    <1+A>   -WN1/2W2 

+  |  N1/2  W]    (1   +  a)  +  N3/?  W,   +  ^  N_1/2  W2  (1   +  2A) 


dNl/2 
dt 


-|N1/2W1    (1   ,A)   -2-N1/2W2  (1   +2A)   -t|n1/2W2 
+  I  N3/2  Wl  +  1  N5/2  W2  +  if  N-3/2  W2  ^   +  2A> 


dN 

"dT^  =  "  I  N-l/2  Wl    -  7  N-l/2  W2  -  TO  N-l/2  W2   ^  +  2A> 

+  |  N_3/2  W]    (1   +  A)   +  \  N_5/2  W2   (1   +  2A)   +  ^  N3/2  W2 


-^  -  -  N.3/2  W]    -  |  N3/2  W,    (1   +  A)   -  ^  N_3/2  W2   (1   +  2A) 

+  N-5/2Wl    <]   +A)+IN-1/2W1   +WN1/2W2 

dN 


'-5/2 


dt 


>5/2  W]    (1   +  A)   -  1  N_5/2  W2  (1   +  2A)  +  N_3/2  W]   +  1  N.1/2W2 
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where  the  downward  transition  probability  per  unit  time  for  Am  =  1  and 
Am  =  2  are  given  by  the  upward  transition  probabilities  per  unit  time 
multiplied  by  the  Boltzmann  factor  1  +  a  and  1  +  2a  respectively,  where 
A  =  ~Yf  and  one  notes  nco.  <<  KT  in  the  temperature  range  of  interest. 

Let  N  5,2,  N  3,2,  and  N°, /2  be  the  number  of  nuclei  in  state 

5   3   1 
2'  2'  2' 

first  approximation  one  would  have 

N 

H/2 

N 
where  nQ  =  g,  and  N  is  the  total  number  of  nuclei.  Similarly. 


m  =  t,  o",  p  respectively,  at  equilibrium  with  lattice  bath.  To  the 


NP  K/2  -  N3/2>  ~-%'~-   N5/2  "  N3/2  ^ 


^  («5/2  "  N?/2>  ~~   2l1o  A  *  N?/2  "  N?/2  ^ 

^1/2 

It  follows  that 

-dp=-  (N5/2  -  N°/2)  W]  -  1  (N5/2  -  N°/2)  W2 

+  (N3/2  -  N°/2)  W1  +  1  (N1/2  -  N°/2)  W2  .  (A-28) 

Defining  nm  =  N  -  N^  ,  one  can  write  the  following  equations  from  (A-25): 

~dF  =  -n5/2  Wl  "  \  n5/2  W2  +  n3/2  Wl  +  \  nl/2  W2  ^"29) 

3/2     2  9 

"St-  =  "  5"  n3/2  Wl  -  n3/2  Wl  "  TO-  n3/2  W2  +  M5/2  Wl 

+  I  nl/2  Wl  +  if  nl/2  W2 
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dF-  =  "  I  nl/2  Wl  "  1  nl/2  U2  "  ll  nl/2  W2  +  I  n3/2  Wl 

1  9 

+  2   n5/2  W2  +  TO  n-3/2  W2 
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+  2"rl-5/2W2     +Wn3/2W2 
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"dt" 
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By  multiplying  each  equation  with  the  corresponding  m  and  summing  all 
togather,  ens  obtains  the  following  equation: 
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lo  a  first  approximation  and  assuming  that  the  spin  temperature  can  be 
defined,  one  notes 
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ha.  •  ° 

where  A'  =  77^-   and  T  is  the  spin  temperature  which  is  generally  greater 
KTS       s 


:han  the  lattice  temperature.  One  then  concludes  that 


<Iz>  -  <l2> 


dt 
which  implies 


■  ■  I  ^v  -  ^J  ^1 +  V 


(A- 32) 


1  =  §  (W,  +  w2) 


(A-33) 


and  finally,  it  follows  from  (A-23) ,  (A-23  )  and  (A-33)  that 


r  a2 


T]  "  125 


e2Qq  1 
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in  the  fast  motion  regime. 


APPENDIX  B 

NMR  SPIN-LATTICE  RELAXATION 

THROUGH  THE  QUADRUPOLE  INTERACTION  IN  SOLID  B-CO 


In  solid  e-CO,  the  time  averaged  quadrupcle  interaction  represents 
a  small  perturbation  term  to  the  total  Hamiltonian  and  the  Zeeman  energy 
levels  are  unequal.  In  the  absence  of  spin-exchange  transitions,  the 
Boltzmann  distribution  of  the  populations  among  the  energy  levels  cannot 
be  maintained.  In  such  a  case  the  spin-lattice  relaxation  is,  in  general, 
not  characterized  by  a  single  relaxation  rate.  In  addition,  the  relative 
importance  of  the  different  relaxation  rates  depends  on  the  nature  of  the 
relaxation  processes  and  on  the  initial  conditions  [38].  However,  since 
the  quadrupole  splittings  are  small  compared  with  the  Zeeman  energy  levels, 
this  simplifies  the  problem.  Considering  for  spin  I  =  j,   from  (A-25) 
one  has 

(B-l) 

dt =  <N5/2  "  N3/2)('2W1  "  K>  +  f  (N3/2  "  NI/2>(W1  +  V 


d(N5/2  -  N3/2> 


+  TO  <Nl/2  "  N-l/2>  W2  +  I  V   ^1    "W2: 


d  ^dt"  "V-      =   (N5/2  "  N3/2>    <W1   "  \V~  VS/2  ~  N1/2)(4W1  +  *%> 
+  4  (N-l/2  "  N-3/2}  W2  "  I  noA   (W1   -5W2} 
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'^dt        =  1  (N5/2  "  N3/2>  U2  +  I  <N3/2  "  Nl/2}  ^  "  W2> 

-f  (N1/2-  N_1/2)  W2  +  f  (N.1/2  -  N_3/2)(W1  -W2) 

d<N-l/2  -   N-3/2}   =   <N-3/2  "  N-5/2><Wl    "  1W2} 
dt 

"  I  <N-l/2  ~  N-3/2>   <4W1  +  7W2) 

+  TO  <N3/2  "  Nl/2>  W2  "  I  V  {W1   "  5W2) 

d(N  ^/9  -  N  -/9)  , 

'3/dt        '      "  <N-3/2  "  N-B/a"'2"!   "  ?W2> 

+  I  <N-l/2  "   N-3/2>    <W1   +  V 

+  TO  ("1/2  "  N-V2>  W2  +  5  V   (W1   -  H2> 

where  W, ,  W„  and  A  have  the  same  definitions  as  in  Appendix  A.  If  one 
now  defines  N  .,  ,p  =  N^i  -  N    from  the  symmetrical  equations  above,  it 
follows  (B-2) 

N2  =  -N2  (2W1  +  1  W2)  +  |  N]  (^  +  W2)  +  ^  NQW2  +  |  nQA  (W]  -  J  W2) 

N1  =  N2(W1  -  g-Wg)  -  \  N1(4W]  +  7W2)  +  ^  N_1Wg  -  J  nQA  (W]  -  5Wg) 

NQ  =  1  N2W2  +  |  N]  (W,  -W2)  -  f  NQW2  +  f  N.,  (W1  -W2) 

+  I  N  ,W,  +  4  n  A  (-W,  +  2W9) 
2-2  2   5  o     1     2 

with  two  symmetric  equations  for  N  -,  and  N  ?. 
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The  solutions  of  these  three  equations  are  characterized  by  only 
three  relaxation  times.  To  simplify  the  equations  in  (B-2),  one  defines 

Np  =  Np  -  nQ  (B-3) 

where  n  and  a  have  the  same  meanings  as  in  Appendix  A.  (B-2)  then 
becomes 

N2  =  -NJ  (2W]  +  1  W2)  +  |  Nf  (W1  +  W2)  +  ^  N5  W2  (B-4) 


N1  =  N^  (W]  -  ^W2)  -  \-  Nf  (4W1  +  7W2)  +  ^  N^  W2 

NQ  =  1  N2W2  +  f  N;  (W1  -  W2)  -  f  N6  W2  +  f  N:i  (W,  -  W2)  +  1  N2W2  . 

Now  assuming  a  precessional  motion  in  3-CO  and  for  simplicity  an 


axially  symmetric  field  gradient,  one  can  calculate  |Vrj2  from  (A-17), 
(A-17'),  (A-13)  and  (A-18')  for  i  =  ±1  and  i  =  ±2 


9 


V^l2  =  J  sin2  b  cos2  p  V2Z  (B-5) 


|v;2!2  =resin"  3  vzz  ■ 

In  chapter  III  it  is  shown  that 

3  cos2  8-1*0 
so  that,  (B-5)  can  be  reduced  to 


iv±l  I2 

=  1V2 

2  VZZ 

W[\* 

=  lv2 

4  VZZ 

1*6   i2 

=  0 

(B-6) 


where  the  V "tennis  are  the  EFG  tensor  components  in  the  crystal  frame, 
while  the  tensor  transformation  from  the  crystal  reference  frame  to  the 
laboratory  frame  follows  the  same  procedure  as  in  (A-17)  and  (A-18). 


V4l  =  V2  D21  +  V:2  D-21  +  Vl  Dll  +  V-'l  D-ll 

V+2  -  T2   D22  +  V;2  D_22  +  Vf  D12  +  V:i  D.i2 


(B-7) 
(B-7') 


For  a  powder  sample  using  (A-16)  and  (A-18)  and  taking  over  all  the 
possible  directions,  one  calculates 


V2 

v77 


H  !    10  "Z 

V 


IV   |2  -  _1  V2 
,v+2'    10  ;ZZ 


(B-8) 

(B-8') 


The  same  equations  will  hold  for  Iv_-j|2  and  [V_2I2  ,  assuming  again 


Vi  (t)  V|  (t  -  t)  =  |Vi  (t)!2  e 


for  i  =  ±1 ,  ±2  where  xX  characterizes  the  correlational  precessional 
motions.  Using  (A-6)  and  (A-6')  one  obtains 

T2 


u  =  \*2M       48 
1   ^  40TiJ 


JL  T2 


W„ 


fe!Qgl2  48 ~2— 


(B-9) 
(B-9') 


From  (B-9)  and  (B-9')  it  is  obvious  that  W,  =  W«.  In  the  fast  preces- 
sional motion  regime  u,  xX  «1 .  It  follows  from  (B-9)  and  (B-9')  that 


W]  =  48 


W2  =  48 


re2Qg^ 


h40, 

^e 2_Qol  2 

ti40 


(B-10) 
(B-10') 


That  is,  w\  =  W?  in  the  fast  motion  regime. 
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Now  from  (B-4)    it  is  possible  to  write 

N'     =  _  f.  U  N-  +  —  W  N'  +  -- -  UN'  (B-ll) 

,N2  2   T7       5  Tl       10     10  K         ' 

^0        2    TN2       5     1H0       2     1-2 

Since  one  considers  the  initial  condition  that  all  lines  are  saturated, 
it  follows  from  the  symmetric  equations  that  N'-,  =  NC  and  N'2  =  N«. 
(B-l'l)  is  further  reduced  to 

N2  =  "  I  W1N2  +  f  W1N1  +  TO"  W1N6  (B"12) 

W'  =  1  W  N'  -  —  W  N' 

"i     2    iIN2     io   ri 

Q 

N'   =  W  N'   -  —  W  N* 

M0        1N2      5  wro 
The  general  solutions  of  (B-12)  will   be 

where  p  =  0,  1,  2.     Substituting  (B-13)   into   (B-12)  and  comparing  the 

-Ait         -  A  o  t  i      - X a t   j.  i 

e     x    ,  e     *■     and  e     d     terms  one  has 

(-  §w1  +  \x)  an  +  |Wj   an  +  ^  W]  a01     =     0  (B-14) 

^■w1  a21  +  (-  ||  W,  +  Ax)  ail     =     0 

W]  a2i  +  (-  |w]  +  A!)  a0i     =    0 

and  similar  equations  for  e~  2°  and  e~   3  terms. 
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In  order  that  all  a-,  f   0  for  i  =  0,  1,  2;  j  =  1,  2,  3  one  needs 

1  %j 


-  2  Wl   +  X 


JrW, 


5     1 


-iW 

10  wl 


"Wwi+A  ° 


0  -  l-W,  +  A 

o     1 


:b-i5) 


=  0 


This  is  reduced  to 


(5A  -  4W.)(2A   -  3Wj)(lQX  -  33  Wj)   -  0 
i.e. ,  with  solutions 


4  3  33 

X  =  5  Wl   •     2  Wl    '  TO  Wl 


B-ie: 


In  the  case  that  all  the  lines  are  saturated  at  time  zero,  one  has 
from  (B-3) 


and 


N2  -  N1  -  NQ  -  0 


N^  =  N,  =  Ni  =  -n  A 
c         1    0    o 


(B-17) 


4         3  33 

Let  \1  =  -F- U-j  ,  x2  =  o"  V7-,  and  ;\3  =  tq  W-, .  Comparing  with  terms  of 

e"At,  (B-12)  gives,  at  t=  0 
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(B-18) 

(-  ~  a2i  -  a22  +  jf  a23)  +  |  (an  +  a12  +  a13)  +  -^  (ao:  +  a02  +  a03)  =  0 

i  i  i 

2  (32i  +  322  +  a23/  +  (-  J  an  +  5  ai2  +  2a13)  =  0 

3      3 

(a2i  +  a22  +  a23)  +  (-aoi  -  yg-  ao2  +  jf  aos)  =  0. 

From  (B-17)  and  (B-13)  one  has  at  t  =  0 

-  n  A  =  a2i  +  a22  +  a23  (B-19) 

-  nQA  =  an  +  ai2  +  ai3 

-  nQA  =  a0i  +  ao2  +  a0  3 
(B-18)  and  (B-19)  enable  one  to  deduce 

a2i  =  a0i  =  an  (B-20) 

3  2 

a22  "  TO  30  2  "  "  5  a12 

a„3  =  _4  a,3  =  |  a23 
which  gives 


a21  -   aIl  _  a0i  "  _nQA 

and  all  other  a-.  =  0  for  i  =  0,  1 ,  2  ;  and  j  =  2,  3.  The  solutions  will 
then  be 

Np  =  nQA  (1  -  e  -V5XW,tj  (B.22) 


for  p  =  0,  1,  2. 


Therefore,  in  the  case  of  complete  saturation  and  W,  =  VL,  the  spin- 
lattice  relaxation  is  governed  by  a  single  exponential  decay  having 


1  -  4  TT 
T,  "  5W1 


(B-23) 


From  (B-10)  and  (B-23)  one  has  the  same  equation  as  in  the  liquid  phase 
in  (A- 34) 


1  _  3 
T1   125 


r  e£Sa 


;B-24) 


where  x'  is  the  correlation  time  for  the  molecular  precessional  motions 
around  the  crystal  C  axis. 


APPENDIX  C 
L1NEWIDTHS  AND  MOLECULAR  MOTIONS  IN  a-CO 


One  assumes  that  the  molecules  in  a-CO  execute  librational  motions 
and  sudden  reorientations.  The  tensor  components  in  the  laboratory  frame 
and  those  in  the  principal  reference  frame  are  related  by  the  rotational 
operator  R  (a,  <f>x) 


]/..  =   R"1  (a,  k)  V(j  R  (a,  ♦*) 


(C-l) 


where  R  (a,  <f>.)  =  R  ($.)  R  (a).  The  Vr.  and  V..  terms  are  the  EFG  ten- 
y     x         ij     tj 

sor  components  in  the  principal  and  laboratory  reference  frames,  respec- 
tively. $1  is  the  out-of -plane  librational  angle  end  a  is  the  reorien- 
tations! angle. 


V{.  =  eq 


£00 

o  - 1   0 

0    0    1 


;c-2) 


In  (C-2)  a  symmetric  field  gradient  has  been  assumed.     It  so  happens 
that  by  assuming  a  small  angle  for  <$>±i 


v„  -  -  2§  0  -  3*5  ) 

V      =  -  -|  (1  +  3+1  sin2  a  -  3  sin2  a) 

Vzz  =  "  ~2  (1   +  3<f>?t  cos2  a  "  3  cos2  a) 
Vxy  =  "  ~2     3K  sin  a 


(C-3) 
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89 


v.,  ■  sa   3K 


X2 


COS  a 


Vw,  =  --§■  3  sin  a  COS  a 
yz    ^ 

For  a  sudden  reorientation  of  an  angle  from  0  to  u  as  is  the  case 
discussed  in  NaNCL  by  Ambrosetti  et  al .  [57]. 

«■  (f  +  *£)  -  (f  +  *£)  s(t)  +  ♦„  (C-4) 

where  s(t)  may  take  values  +1  or  -1,  and  <f>*  is  introduced  to  account  for 
jumps  not  exactly  equal  to  it.  <j>„  is  the  in-plane  librational  angle.  By 
assuming  small  angles  for  <j>„  and  <J>*  and  expanding  sin  a  and  cos  u,  one  has 


*£)<-$*5.(t»  +  *.  [*£+<i -$s(tn  l    ' 


K*2 


sin2c^  (1  -  s(t))  V+  ♦«  +  4n£(l   -  s(t)) 

COSa^    (1    -*£    [*J+(1    -*£)    S(t)j+M-  (1    -S(t)) 


C0S2a   -    1    +  V  (S(t)    -    1)    -  *4-  (1    -    S(t))    -    6„2 


From  (C-3)  and  (C-5)   it  follows  that 

vxx  -  -  e-f  (1  -  3*|)  (c.6) 

vyy  ■  -  *£  [i  -  ^r  (]  -  s^))  -  ™> +  -  Mr  (1  "  s(t))] 

V__  =  .  Q  [-2  +  3*5  +  3*2  +  3iji  (1   _  s(t))  +  3?^  (1   _  s(t))] 


\y   =    -^**     (-^  +  ^S(t)    +   *„    S(t)) 

vxz=      ^U^+S(t)(l   -*Ji)+iif  (1   _s(t))] 


v.-«-^(-*Js(t)  +*£+♦;,)  . 


yz 
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Assuming  i>H  and  $*  are  not  correlated  and  also  that  <  $„  >  =  0,  the 
librational     average  of  (C-6)   is  given  by 

<  Vxx  >  =  -  3J.  (T-  3  <  9,2  >  )  (C-7) 

<  V       >  =  -  S|  [1   -  3^i  (1   -  <  s(t)   >)   -  3  <  9„2  >] 

<  V       >  =  -  £§.  [_2  +  3<^2>+  3  <  <t»„2>+  ^  (1   -  <  s(t)   >)] 


<v     >=<v       >=<v     >=0 
xy  xz  yz 

The  quadruoole  Ha/nil  Ionian  can  be  divided  into  two  parts 

HQ  =  *  HQ  >  +  [AQ  "  *  HQ  >]         *  (C"8) 

The  average  Hamiltonian  <  hq  >   gives  the  NQR  frequency,  while  the  second 
part  drives  the  relaxation  process  and  causes  the  line  to  broaden  [58]. 

The  fluctuating  terms  of  the  components  of  the  EFG  tensor  can  be 
written  as 

v7Z  -  <  vxx  >  ■-    2fa  ($x2  .  <  9x2  >  )  (c-9) 

vyy  -  <  vyy  >  =  -3f  [♦;  -  <  9„2  >  -  ^f?-  (s(t)  -  <  s(t)  >) 

-^    (1  -  s(t))] 

VZ7     -     <     V^     >     =     -    ^    |>A2     -     <     K2    >+     9n2     -     <     9„2     > 

+  i£(s(t)  -  <  s(t)  >  )  +2*j£(l  -  s(t))] 

Vxy  "  *  Vxy  *  =  "  ^  ♦*   C  "  *?  (1   "  S(t))  +  ♦"   S(t)] 
Vxz-<Vxz>a      ^♦xs(t) 

v    -  <  v    >  -  -  3§a[4o  -  s(t))  +  *^ 

yz  yz  c         c 
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where  in  (C-9)  consideration  is  given  only  to  the  second  order  of  $„, 
<j>j,  and  $*. 

For  spin  I  =  j   ,  the  NQR  frequency  between  m  =  ±j  and  m  =  ±j   ,  T2 
is  given  by  [58]. 


1.1+11  +  1 

T2       T2       2      I  T±5/2       T±3/2  J 


(C-10) 


where 


r'(t)  u*(t  -t)  dx 


ii'(t)  is  the  departure  of  the  instantaneous  resonance  frequency  from  its 
average  value  and 


5/2 


T±3/2   m*  ±3/2 


itf  ±5/2   W±5/2'  m 
I  w  3/2,  m 


(Oil) 
(C-12) 


where  T+5,2  and  T+3/~  are  the  lifetimes  of  the  states  m  %c/?  an'J  m  %o/? 

respectively.  From  (A-l )  and  (C-9)  one  will  have,  assuming  order  phase 
<  s(t)  >=  0 


<°'(t) 


9e  Qq 
2oti 


/6 


[*i2  -  <  <$>±2  >  +  <J>u2  -<  4>h2  >         (C-13) 


i*2 


+  ^-    s(t)   +  I*,,**  s(t)] 


Since  $„  and  ^  occur  at  higher  frequencies,  for  simplicity  one  may  neglect 
these  terms.  One  would  then  have 


r(t)o>-(t  -t)  dr  = 


9e2Qq 

2dfi 


2    ,   f. 


s(t)  s(t  -t)  di  (C-14) 


Jo 
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Defining 


s(t)   sTT^T]"    e~1ClJT  dr     =     <]„    (w) 


it  follows  from  (C-14)  and   (C-15)  that 

f     "l2 
1  .   9e  Qq 

T2   [  20n 
Consider  now  for  -■         terms.  From  (C-ll),  one  writes 


3  <  <j>*k  > 
16 


Jn  (0) 


±5/J 


+  w±^/: 


T±5/2        ±5/2>±3/2        ±5/2, ±1/2 
Using  (A-l),   (A-4),   (A-4'h    (A-6)  and  (A-6')  one  obtains 


f±5/2       (  40"n  t 


2     r 

80    I 

L 


J-o 


T^tTlt^t  -t)  e"1WT  dx 


_1Q 

2 

40 

1  40n    J 

L_           . 

V^(t)  V*2(t  -TT    e 


-  1  ?WT 


f. 


+  V  9(t)   V*  (t  -tT    e 


■1*MT 


-2V 


di 


Now  one  needs  to  consider 


where 


A  = 


r 

CO 

A  •--• 

r    ■ 

3eq_ 
"2 

J-oo 

f 

2  r 

LI. 

lU)T 


*2~ 


V±1(t)  V^  (t  -T)     e",a)T  dx 

00 

CTt)   s(t)  k  (t  -t)   s   (t  -  tT  e 

-00  . 


-lOJT 


+  |     *J-s(t)  s(t  -t)  e-1UT  dT  + 


(C-15) 


(C-16) 


V+1(t)   V^Tt~^T  e"ia'T     dr  (C-17) 


(C-18) 


-T(OT 


$„(th»[t  -  x)  e   U0T  dr 


J. 


and  where  no  correlation  between  the  value  of  <f>*  and  s(t)   has  been 
assumed. 


«' 


If  one  assumes  an  harmonic  oscillator  for  ^  and  6.,,  and  for  s(t) 
:ponential   cc 
from  (C-iS)   that 


an  exponential   correlation  function  with  characteristic  time  xX\  it  follows 


A  = 


2§a 

2 


2       2tX' 
<  4>x     >       ^ 


(C-19) 


i*2 


1   +  (wr~+  oj)2t  — : 
2 


r+  A-2  J 


One  notes  that  o^  >>  u    and  iaxrX^»     1 


I       J 


t,*2 


"1     +    a)"2""^ 


(C-20)  can  be  reduced  to 


A  -   f3e3.        <  ft*2 
\  "  J 


TT~3^f 


:c-2D 


While  for 


v;9TFj  V*Tt  ~T7    e"1^1  dx 


±z 


Since  V^    involves  V+«»  one  needs  only  to  consider,  assuming  <s(t)>  =  0, 

vyy  "  <  v'yy  >  '  "  ~f  [  -^r1—  s!t»  "  %^  i1  "  «<*»      (c-22> 

and 


V       -  <  V 
xy  xy 


[■ 


s(t))  +  K4-H  s(t)] 


If  one  further  assumes  that  no  correlation  exists  between  <+>„,  ^,  (j>* 
and  s(t),(C-22)  can  be  simplified  as 


V   -  <  v   >  .  -  M  3<|>*2   sft) 

yy         yy  z  2     su; 


(C-22-) 


V   -  <  V   >  «  0 
xy     xy 


35 


It  follows  that 


3  - 


9  <  <b*k  > 


2   i 


iTtTTr^Tye"1^  dx  (C-23) 


3eq 


2     9  <  **h  >       .       {  3     x 
— 4 JH    (  2  u) 


From  (C- 17) ,   (C-21)  and  (C-23) 


•±5/2 


{     80n    J     i_  2 


(C-24) 


1 


Considering  now  the  y  term  and  using  (A-l),   (A-4)»   (A-4')»   (A-5) 

'±3/2 
(A-5'),    (A-6)  and  (A-6'),  one  writes 


3/2 


40h 


L 


80 


V+1(t)   V^(t  -  t)  e",aT  dT 

r^rv^Tt'^T  e-iuT  dx  1 


(C-25) 


4dti 


V+1(t)  V*^Tr^~TTe"'2T  dx 

.0) 

V_1Ttrvf1Tc  -  xT    e"j2T     dx 


h   f  eO   |2 


40n  J 
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UT       . 


V+2(t)  V*2(t  -  t)  e"'2T  d 


-1*T 


V_2(tTV*2(t  -  x)     e-1?1     dx 


Employing  the  A  and  B  terms  as  discussed  above,  one  writes 


±3/2 


9e2Qg 
20n 


+  f  <    **2    > 


T2' 


<  *       >  t~+ — z* — ^T 


1   +  jf  w2  x£ 


-    +f<^>  Jn(  1.) 


(C-2i 


Rfi 


From  (C-10),   (C-24)  and  (C-26),  one  has 
3 


1 

V 

f  9e2Qg_  I2 
I     2<A     J 

^ 

-  !  <  **" 

i 

-  1  <  **2 

tI <  ♦** >  je (o)  +  £<♦**>  Jn(j«) 


1    (3  \      4.  10 


1  +  1  CO2  t- 


,*2 


1  +  U)Z  T^'T 


(C-27) 


If  one  assumes  fast  motion,  i.e.  Jfl  (0)  -  Jj  (i  u)  =  J-  (w)  *  j"n  (|  us), 
one  can  neglect  terms  involving  <  <j>*'+  >  since  <t>*  is  a  small  angle.  One 
then  has 


(C-28) 


9e2Qc 

20h" 


12 


s*2 


T2' 


TO 

9   1  +  (^  T^ 


x2 


^7: 

4  T2 


These  terms  on  the  right  hand  side  of  (C-28)  came  from  the  lifetime  of 

5        3  -i 

trie  states  m  =±-^-and  rn  =  ±j,   while  the  adiabatic  term  4.  has  been  neg- 
lected. One  can  also  show  that  for 


r£*  -  1.1! 


(C-29) 


(C-28)  has  a  maximum. 

One  thing  worth  noting  in  (C-7)  is  that  <  V   >  shows  the  usual 
Bayer-type  dependence  for  the  NQR  resonance  frequency  even  though  it  is 
derived  from  the  small  angle  $*.  One  also  notes  that  for  large  $*,  i.e. 
significantly  different  from  tt  reorientation,  <  V   >  will  be  affected. 

In  the  fast  motion  regime,  u)t~  «1 ,  (C-28)  can  be  further  reduced  to 
1  =  4  f  9e2Qq  ) 2 


20ti 


<   <j>*"    >      tX- 


(C-30) 


which  relates  the  linewidth  and  the  correlation  time  of  the  molecular 
reorientational  motion. 
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